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Expressions are given for the diffraction field far behind, and the surface 
currents on, a parabolic cylinder. Approximate values for the field strength 
and current density are given when the radius of curvature of the cylinder is 
large compared to a wavelength. The fornmlas may have value in predicting 
the shadows that are cast by hills in microwave propagation. The idea of 
representing hills by knife-edges has been used successfully by a number 
of investigators. The theory of the parabolic cylinder indicates that such a 
representation is valid even for gently rounded hills when the angle of diffrac- 
tion is small. On the other hand, when the angle of diffraction is so large 
that the knife-edge calculations do not apply, the results presented here may 
be used. 

1. INTRODUCTION 

A number of investigators have studied the effect of hills on the 
propagation of short radio waves. Experiment has shown that the field 
far behind a hill may be computed, to a reasonable degree of accuracy, 
by assuming that the hill acts like a knife-edge (half-plane) 1 . The ques- 
tion naturally arises as to the conditions under which such an assump- 
tion is permissible. Here we attempt to throw some light on this ques- 
tion by taking the hill to be a parabolic cylinder. 

Our results indicate that, for small angles of diffraction, even gently 
curved hills act like knife-edges. However, for larger angles correspond- 
ing to points deep in the shadow or to points high in the illuminated re- 
gion, it may be necessary to use the more exact formulas which take the 
curvature of the hill into account. 



1 See, for example, Ultra-Short-Wave Propagation, J. C. Schelleng, C. R. 
Burrows and E. B. Ferrell, Proc. I.R.E., 21, pp. 423-463, 1933. 
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As an application of our results, a brief study is made of the hypo- 
thesis that very short radio waves are transmitted far beyond the hori- 
zon by successive diffractions over hills or ridges. The ridges are as- 
sumed to be of equal height, to be 40 miles apart, and to have a radius 
of curvature of 100 meters at their crests. At the frequencies considered 
(30 and 300 mc) and at the small angles of diffraction required to go 
from one crest to the next, the ridges act like knife-edges. 

At 30 mc and a distance of 280 miles the calculated field is in fair 
agreement with the observed field 2 . At 300 mc and at the same distance 
the calculated field is about 50 db below the observed field. This sug- 
gests that the 30 mc long distance propagation may possibly be ex- 
plained by successive diffractions. The discrepancy at 300 mc may be 
due to any one of a number of reasons. For example, it may be due to 
the effect of the non-uniformity of the atmosphere, or to the roughness 
of our approximations (for one thing, we neglect reflections from the 
ground between the ridges). 

The first theoretical work on the diffraction of plane electromagnetic 
waves by a parabolic cylinder apparently was done by P. S. Epstein. 3 
His work makes use of a series of parabolic cylinder functions. When the 
cylinder is large many terms are required for computation. By "large" 
we mean that the radius of curvature at the vertex of the cylinder is 
large compared to the wavelength of the radio wave. 

An entirely different approach was used by V. Fock 4,5 in 1946. In the 
first paper Fock sketches the derivation of an integral for the current 
density on a large paraboloid of revolution. In the second paper he re- 
derives this integral by considering the form assumed by the field equa- 
tions when a plane wave strikes a gently curved conductor at grazing 
incidence. His result gives the change in current density on a large and 
highly conducting parabolic cylinder as we go from the illuminated 
region into the shadow. 

In 1950 K. Artmann 6 examined the diffraction field far behind a large 
circular cylinder. He showed that, for small angles of diffraction, the 

2 A summary of experimental data is given by K. Bullington, Radio Trans- 
mission bevond the Horizon in the 40-400 Megacycle Band, Proc. I.R.E., 41, pp. 
132-135, 1953. 

3 Dissertation, Munich (1914). A more accessible account of this work is given 
in the Encyklopaedie der Math. Wiss. 5, Part 3 (1909-1926) Phys. p. 511. Sec also 
II. Bateman, Partial Differential Equations of Math. Phys., (Cambridge Univ. 
Press 1932) p. 48S. 

* The Distribution of Currents Induced by a Plane Wave on the Surface of a 
Conductor, J. Phys. (U.S.S.R.), 10, pp. 130-136, 1946. 

6 The Field of a Plane Wave Near the Surface of a Conducting Body, J. Phys. 
(U.S.S.R.), 10, pp. 399-409, 1946. 

6 Beugung polarisierten Lichtes an Blenden endlicher Dicke im Gebiet der 
Schattengrenze, Zeitschr. fur Phys. 127, pp. 468-494, 1950. 
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diffraction pattern is shifted by an amount proportional to the J-^ power 
of the radius of curvature. Whether the shift is towards the shadow or 
in the opposite direction depends upon the polarization of the incident 
wave. 

In this paper we derive some of the results given by Fock and Art- 
mann by starting with Epstein's series. In addition we investigate the 
diffracted field at a great distance behind the cylinder. The cylinder is 
assumed to he a perfect conductor in all of our work except for a few 
equations given near the ends of Sections 4, 6, and 7. The procedure is 
similar to that used in the smooth-earth theory. 7 The series is converted 
into an integral and then the path of integration is deformed so as to 
gain as much simplification as possible. As might be expected, the results 
for a large parabolic cylinder are similar in some respects to those for a 
smooth earth. Much of the work requires a knowledge of the behavior 
of parabolic cylinder functions of large complex order. Although several 
studies of this behavior have been published, the results are not in the 
form required. For this reason, and for the sake of completeness, several 
sections dealing with the properties of parabolic cylinder functions have 
been included. 

Incidentally, W. Magnus 8 has studied the field produced by a line 
source located at the focus of a parabolic cylinder. However, his problem 
is somewhat different from the one with which Ave are concerned. 

I am grateful to Prof. Erdelyi of the California Institute of Tech- 
nology and to my colleagues at Bell Telephone Laboratories for helpful 
discussions and references which resulted in a number of improvements 
throughout the paper. I am also indebted to Miss Marian Darville for 
performing the rather laborious computations upon which the various 
curves and tables are based. 

2. DISCUSSION" OF RESULTS 

Various expressions are given later for the electromagnetic field in 
terms of parabolic cylinder functions. In this section we shall confine a 
good share of our attention to the case in which the cylinder is very 
large compared to a wavelength so that the cylinder functions may be 
approximated by Airy integrals. As in the remainder of the paper, we 
shall be concerned chiefly with the field behind the cylinder and the 
current density on the cylinder. 

7 By "Smooth-earth theory" we mean the formulas for the field produced by a 
dipole near a large sphere. A complete discussion of the theory is given in the 
book by II. Bremmer, Terrestrial Radio Waves (Elsevier, 1949). 

8 Zur Theoric des zvlindrisch-parabolischen Spiegels, Zeitschr. fttr Physik, 
118, pp. 343-356, 1941. 
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It turns out that the results for the parabolic cylinder are closely 
related to those obtained by Sommerfeld 9 for the diffraction of a plane 
wave by a perfectly conducting half-plane. In fact, the two fields are 
surprisingly similar in the region of the shadow boundary. More pre- 
cisely, the fields are similar for values of the angle \f/, defined in Fig. 2.3, 
such that (roughly) 



\f/ in radians | < - 



wavelength 



i /a 



radius of curvature of cylinder at its crest_ 

where the coefficient H fias been selected somewhat arbitrarily. For 
larger values of | ^ | the difference between the fields becomes pro- 
nounced. As we go deeper and deeper into the shadow, i.e. as \f/ becomes 
more and more negative, the field behind a cylinder ultimately decreases 
exponentially with \f/. On the other hand, the field behind a half-plane de- 
creases roughly as 1/| ^ |. Since the exponential function decreases more 
rapidly than does 1/| \f/ |, the shadow behind a hill is darker than the one 
behind a half-plane. High in the illuminated region the field consists of 
the incident wave plus the wave reflected from the illuminated portion 
of the cylinder. For the half -plane this reflected wave is negligibly small 
until \p reaches 180°. 

First we shall review the situation pictured in Fig. 2.1. An incident 
wave comes in from the left and strikes a perfectly conducting vertical 
half-plane which casts a shadow as shown. The electric and magnetic 
intensities are proportional to exp {iwt) where t is the time and w is 
the radian frequency. The unit of length is chosen so that X, the wave- 
length, is equal to 2ir. This is done in order to simplify the expressions 
we have to deal with. For example, a plane wave of unit amplitude 
traveling in the positive x direction, as shown in Fig. 2.1, is represented 
by exp (— ix). 

Sommerfeld's exact expressions, for the special case of horizontal 
incidence shown in Fig. 2.1, may be written as 

(hp) E = (e- ix + S 1 ) + £ 2 (0), (2.1) 

(vp) H = (e~" + Si) + £ 3 (0), (2.2) 

where (2.1) holds when the electric intensity E is parallel to the edge, 
and (2.2) when the magnetic intensity H is parallel to the edge. From 



9 Math. Annalen, 47, p. 317, 1896. Sommerfeld's results have been described in 
a number of texts on optics. The book, Huygcns' Principle by Baker and Copson 
(Oxford 2nd edition, 1950) deals with this and many similar problems. See also 
Chap. 20 of Frank-von Mises, Die Differenlial-und Inlegralgleichungen der Me- 
chanik und Physik, 2nd edition, Braunschweig: F. Vieweg and Sohn (1935). 
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the analogy with the radio case, these two polarizations will be termed 
"horizontal polarization" (hp) and "vertical polarization" (vp), re- 
spectively. The incident plane wave for hp is assumed to have an E of 
unit amplitude. This is indicated by the exp (-ix) in (2.1). For vp the 
incident wave is assumed to have an H of unit amplitude. The S's 
are denned by the Fresnel integrals 



c- ix + Si - (*/*)' 



L" 

J— 00 



<r" dt, 



flk- -w*) w *~* /V'" dt, (23) 

&(0) = -5,(0) - -(t/*) w *+* /"•"*"*, 



fa = (2r) 1/2 sin|, 



f, = (2r) 1/2 cos£, 



i 1 ' 2 = exp (tV/4), 



where (r, tp) are the polar coordinates shown in Fig. 2.1 [S\ and £ 2 (0) 
= - S 3 (0) for an arbitrary angle of incidence are given by Equations 
(5.3), (5.6), and (5.20) of Section 5]. We use the notation <S 2 (0), »S ;t (0) to 
indicate that these functions are special cases of the functions 82(h), 
Si(h) which appear in the analysis for the parabolic cylinder. 

The principal part of the field far to the right of the half-plane, where 
x is positive, is given by exp (—ix) + Si whose absolute value is plotted 
in Fig. 2.2. The function Si almost cancels the incident wave in the 
shadow, and then drops down to small values outside the shadow. 
The function >S' 2 (0) is always small in the region we shall consider. It 
becomes large only when <p exceeds -k. It then corresponds to the wave 
reflected from the front (left-hand side) of the half-plane. 

When we are far enough away from the shadow boundary to make 



INCIDENT WAVE 



TRACE OF_, 
HALF-PLANE 




Fig. 2.1 — Diffraction of a plane wave by a half-plane. 
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Fig. 2.2 — The approximate value | e~ ix + <Si | of | E j for hp and of \ H \ for 
vp at a great distance r behind a half-plane. Here, as in all of our work, the wave- 
length X is 2ir. For an arbitrary wavelength replace r by 2t//X, etc. 

<p r >>> 1, exp (— ix) + & has the asymptotic expressions [see Equations 
(7.7) and (7.8)] 



& 



c(r)/2 sin - , <p < (shadow) 

a*"* 4- c(r)/2 sin | , <p > 

c(r) = z 3/2 (2 7 rr)- 1/2 e - > . 



(2.4) 



(2.5) 



These expressions lead us to picture the field to the right of the half- 
plane as the sum of the incident wave and a wave, c(r)/ip, spreading out 
from the diffracting edge 10 (for the small p's of interest, 2 sin <p/2 « <p 
even though y>r» 1). In the illuminated region these two waves inter- 
fere with each other to give the oscillations around unity shown in Fig. 



10 See, for example, R. W. Wood, Physical Optics, 3rd edition, p. 220 (Mac- 
.Millan, 1935). Curves of equal phase and amplitude have been computed and 
plotted by W. Braunbeck and G. Laukien, Einzelheiten zur Halbebenen-Beugung, 
Optik, 9, pp. 174^179, 1952. 
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2.2. In the shadow only the edge wave is present and there is no inter- 
ference. Si is not the edge wave in the shadow. 

The edge does not radiate uniformly in all directions. The <p in the 
denominator of c(r)/<p indicates that the edge sends out its strongest 
wave in the direction of the shadow boundary where tp = 0. This ac- 
counts for the decreasing size of the oscillations in Fig. 2.2 as <p becomes 
more and more positive. It likewise accounts for the steady decrease as 
<p becomes more and more negative. 

That jS 2 (0) is small in comparison with exp (— ix) + Si follows from 



S 2 (0)~c(r)/2cos| 



(2.6) 



and the fact that this is small compared to the c(r)/2 sin (<p/2) in 
(2.4) when <p is small. 

So far, we have been discussing a special case of Sommerf eld's results. 
Now we turn to the case of the perfectly conducting parabolic cylinder 
shown in Fig. 2.3. Here, as in Fig. 2.1, the incident plane wave exp (—ix) 
comes in from the left. The fields for the two kinds of polarization are 
given by 



(hp) 
(vp) 



E = (e~ ix + SO + S 2 (h), 
H = (e~ ix + SO + S 3 (h), 



(2.7) 
(2.8) 



where, just as in the half -plane case, the fundamental vectors E and H 
are perpendicular to the plane of Fig. 2.3 and the incident waves are of 
unit amplitude. 

The [exp (—ix) + Si] in (2.7) and (2.8) is exactly the same Fresnel 
integral (2.3) as for the half-plane. S 2 (h) is a rather complicated integral 
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Fig 2 3 — Coordinates used in the discussion of the parabolic cylinder. The 
coordinates such as (0, h) refer to (x, y). The origin of coordinates is at the 
focus of the parabola and h is the height of the vertex or crest, above the origin. 
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(obtained by setting the angle of incidence equal to t/2 in equation 
(5.4)) involving parabolic cylinder functions. When h = the parabolic 
cylinder reduces to a half -plane and Si (h) reduces to the value *S?(0) 
appearing in (2.3). The symbol S 3 (h) represents an integral much like 
S?(h) except that it contains derivatives of the parabolic cylinder func- 
tions. As we might expect, S*(h) and Sz(h) behave in much the same 
way as does S 2 (0). In particular, they are small compared to exp (—ix) 
+ Si at the shadow boundary, and their asymptotic expressions anal- 
ogous to (2.6) hold as <p and \f/ pass through zero. 

S 2 (h) and Sz(h) have been put in a form suitable for computation in 
two cases, (1) when h = 0, which is the half -plane case already discussed, 
and (2) when h and r/h 2 are very large. In the second case it is conven- 
ient to introduce new polar coordinates (p, \f/) with their origin at the 
crest of the parabola as shown in Fig. 2.3. In these coordinates a circular 
cylindrical wave spreading out from the crest is asymptotically pro- 
portional to 

c(p) = i m (2ttp)- 1/2 <T ip . (2.9) 



In Section 8 it is shown that 

e « [*-* + sj, - 4el + c( P )h lls 



-3- 



¥(t) + - exp (*>73), (2.10) 



1/3. 



r = tf'V 

is an approximation which gives the field (for horizontal polarization) in 
the region of the shadow boundary far behind a large cylinder. Our t is 
an approximation to the g used there. Here the subscript p on [exp (—ix) 
-f $i] p means that the quantity within the brackets is to be computed 
as though it corresponded to a half-plane with its edge at the crest of 
the parabolic cylinder, so that p, \p are to be used in (2.3) instead of r, <p. 
Also, 

m ) = 2 ~ i/3 r Ai ^ — ( ri/!>Mr ) — 

Jo Ai(u) — iBi(u) 

... . . (2.H) 

*L 

Jo 



., . , Ai(u) exp (iur) du 
Ai(u) + iBi(u) 



where Ai(u) and Bi(u) are Airy integrals defined by equations (13.12) 
and (13.16), and tabulated in reference. 11 In this paper we find it con- 
venient to use the Airy integrals instead of the related Bessel functions 



" The Airy Integral, Brit. Asso. Math. Tables, Part — Vol. B (Cambridge, 
1946). 
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of order J^. As in (2.3), the fractional powers of i are made precise by 
taking i = exp {iir/2). 

Three kinds of approximations have been made in the derivation of 
(2.10), namety those associated with the assumptions (1) that the angle 
^ is small, (2) that h is large, and (3) that r/h 2 is large. The terms in l/rp 
and 1/r do not cause trouble at \f/ = because their infinities cancel 
each other. 

The counterpart of (2.10) for vertical polarization is obtained from 
(2.10) by replacing E by H and ^(t) by ^(r), where the subscript v 
stands for "vertical"; and ¥„(t) is given by (2.11) when Ai(u) and Bi(u) 
are replaced by their derivatives with respect to u. 

Series for ^(t) 4- 1/t and ^„(t) + 1/t which converge for negative 
(shadow) values of r are given by Equations (7.31) and (7.53), respec- 
tively, with g in place of t. Table 2.1 gives values of ^(t) and ¥,(t) 
which were obtained from the series for t negative, and from numerical 
integration of (2.11) and its analogue for t ^ 0. 

When t is large and positive Equations (7.35) and (7.55) show that 



*(r) 



1/r 
1/r 



(iVr) 1 ' exp (-tY712), 
(itt) 112 exp (t» - it 3 /\2). 



(2.12) 



When r is large and negative the leading terms in (7.31) and (7.53) 
give 

*(t) + 1/t i m 2.03 exp [(2.025 + i 1.169)t], 

# v ( T ) + i/ T ^ _ j 1/3 3.42 exp [(0.882 + i 0.509) r]. 

Now that we have expressions for the field what do they tell us? For 
one thing, they may be used to show that the field near the shadow 





Table 2.1 — Values of *(t) and ^„(t) 




T 


1 * (r) | 


arg. *( r ) 


1 *vM | 


arg. *„(t) 


3 


3.13 


-93.5° 


3.16 


+104.3° 


2 


2.21 


-1.8° 


2.80 


192.4° 


1.5 


1.945 


+21.6 


2.44 


211.3 


1.0 


1.715 


+32.5 


1.985 


219.5 


0.5 


1.486 


34.2 


1.522 


218.6 





1.254 


30.0 


1.089 


210.0 


-0.5 


1.030 


22.9 


0.724 


193.7 


-1.0 


0.823 


15.2 


0.459 


167.8 


-1.5 


0.648 


8.48 


0.317 


130.6 


-2.0 


0.511 


+3.79 


0.281 


92.0 


-3.0 


0.338 


+0.12 


0.288 


45.1 


-4.0 


0.250 


-0.12 


0.264 


22.3 


-5.0 


0.200 


-0.02 


0.221 


9.71 
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boundary is almost the same as for a half-plane. Away from the shadow 
boundary, the field in the shadow can be interpreted as a "crest wave" 
which reduces to the "edge wave" for a half -plane. The crest wave de- 
creases as an exponential function of \f/ in the shadow instead of as l/<p 
in the case of a half-plane. In other words it is much darker behind a 
parabolic cylinder than behind a half-plane — and the larger the cyl- 
inder the darker it is. (A glance at Fig. 2.8 shows that this statement 
must be qualified for vertical polarization by requiring the observer to 
be deep in the shadow.) As Figs. 2.7 and 2.8 show, deep in the illuminated 
region the crest wave behaves like the wave reflected (as computed by 
geometrical optics) from the illuminated portion of the cylinder. 

Now we consider expressions for the surface currents. Let J and J v be 
the densities of the conduction currents which flow on the surface of the 
perfectly conducting parabolic cylinder for the cases of horizontal and 
vertical polarizations, respectively. J is parallel to E and is perpendicular 
to the plane of Fig. 2.3 while «/„ flows in the plane of the figure. ./„ is 
positive when the current flows in the direction of increasing x. In Sec- 
tion 6 it is shown that when h is large, J and «/„ are given approximately 
by 



exp (- ix - tV/3) f" \i ' exp (- i lt6 uy) 



„» r —2/3 _ 

f i ex 

vh lli k L Ai (u) - iBi(u) 

( ■ \ - (2<14) 

exp (— iuy) 



J v 



Ai(u) + iBi{u)_ 
i exp (— ix — ^y/3) f" exp (— i~ uy) 



du, 



f[ 



( L - (2 - l5) 

exp (— iuy) 



Ai'{u) + i Bi'{u)_ 



du. 



These expressions are obtained when the relations (13.17) for Airy 
integrals are used in equations (6.1G) and (6.23). Here f is the intrinsic 
impedance of free space. In the rational MKS system which we use 
fo = 120* ohms. The factor f appears in (2.14) but not in (2.15) because 
we assume the incident wave for vertical polarization (H = 1, E = £ H 
= 1207I-) to be 1207T times stronger than the one for horizontal polariza- 
tion (E = 1). The primes on Ai(u) and Bi(u) denote their derivatives 
with respect to u. The parameter 7 depends upon the coordinate x of 
the point at which the current is being observed: 

7 = x/2h m . (2.16) 
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Equations (2.14) and (2.15) hold only in the region of the crest of 
the cylinder and for large h. Under these conditions x + y 3 /S is very 
nearly equal to the distance along the surface measured from the crest 
of the cylinder. 

An expression equivalent to the one for J v in (2.15) has been derived 
and tabulated by V. Fock. 4 

Series for f «/ and ./„ which converge for positive (shadow) values of 
7 are given by equations (6.17) and (6.24). When y is large and negative 
the application of the method of steepest descent to integrals (6.16) 
and (6.23) leads to 



f J « -(x/h)e- ix [1 + t-/4 T 3 +•••], 
./„ « 2e~ ix [1 - i'/4y 8 + •••], 



(2.17) 



in which x/h = 2yh~ . 

Table 2.2 gives values of h m toJ exp (ix) and «/„ exp (ix). The values 
of J for 7 > were computed from the series, and the ones for 7^0 
were obtained by numerical integration of (2.14). The entries for J v 
were taken from the more extensive table given by Fock. 4 In order to 
express his results in our terms it is necessary to use the fact that the 
radius of curvature at the vertex of the parabola is 2h. A change in the 
sign of i is also necessary because the time enters Fock's work through 
exp (— tut) instead of exp (iut). The values shown were checked for 
7 > by the series and for y ^ by numerical integration of (2.15) 

Fock's table shows that by the time y has reached —2 the value of J v 
exp (ix) has become 1.982 at an angle of +1.45 degrees. This is close to 
the limiting value of 2 predicted at y = — » by (2.17). 

It will be noted that, for large values of h, J is smaller than J v by 



Table 2.2 — Surface Current Densities 






A>''fo/exp (i* + »V/3) 


J v exp (ix + iV/3) 




modulus 


Argument in degrees 


mod. 


Arg. 


-1.0 


2.16 


-25.9 


1.861 


-15.43° 


-0.5 


1.38 


-16.8 


1.682 


4-1.52 


0.0 


0.77 


-30.0 


1.399 


0.00 


0.3 


0.515 


-44.8 


1.197 


-6.06 


0.6 


0.327 


-62.9 


0.991 


-14.23 


1.0 


0.167 


-90.1 


0.738 


-26.63 


1.5 


0.066 


-125.9 


0.488 


-42.57 


2.0 


0.025 -161.6 


0.315 


-57.98 


3.0 


0.0033 -230.7 


0.130 


-87.57 


4.0 


0.00043 -298.0 


0.054 


-116.75 
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the order of h when y is of moderate size. It will be shown later that 
the current density decreases exponentially as one moves into the 
shadow, and that its rate of decrease is related to that of the field as 
shown in Fig. 2.6. 

We now take up the detailed discussion of the expressions for the field 
and the current density. It is convenient to consider the current density 
first. When a plane wave strikes a perfectly conducting plane, the sur- 
face current is proportional to the tangential component of H in the 
incident wave, and flows at right angles to it. In the rational MKS units 
we are using, the surface current density is two times the incident 
tangential H. When we consider the illuminated side of a large parabolic 
cylinder and calculate the current density by doubling the tangential 
component of the incident H we obtain the approximations 

toJtt -2x(4h 2 + x 2 T m e-**, 

(2.18) 
J v tt2e 



—IX 



which hold when x is large and negative. When h is very large but x/2h 
small these formulas agree with the leading terms of (2.17) which were 
obtained from our integrals for the current density. 

Expressions for the current density deep in the shadow may be ob- 
tained from the leading terms of the convergent series by letting y 
become large and positive. The exponential decrease is found to be 

I f J I « 1.43/T 1 ' 3 exp (-1.018a*"*'*), 

| J v | » 1.83 exp (-0.44.t7T 2 ' 3 ), 

where the numbers appearing in these equations are associated with the 
smallest zeros of Ai(u) and Ai'(u), respectively. These formulas for a 
large cylinder are roughly similar to those for propagation over a smooth 
earth. The radius of curvature at the crest of the cylinder is 2h. Setting 
this equal to the radius of the earth gives an exponential rate of decrease 
for J and J v which is the same as that over a smooth earth for the 
two polarizations. Of course, the coefficients multiplying the exponential 
functions are different. This agreement is not surprising since the Airy 
integrals are closely related to the Hankel functions of order J^ used in 
the smooth earth theory. 

The surface current densities as a function of the distance h — y 
below the crest for h = 1000 and f or h = are shown in Fig. 2.4. The 
equation of the cylinder shows that h —y = x 2 /4:h so that, from (2.19), 
foJ and J v decrease in proportion to h~ llz exp [— 2.025/T 1 ' 6 (h — y) 112 ] 
and exp [— .88/T 1 ' 6 (h — y) 1 ' 2 ], respectively, far down in the shadow. 
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Fig. 2.4 — The surface current density is plotted as a function of the vertical 
distance below the crest or edge. The curves for /; = 1000 and h = are obtained 
from Table 2.2 and equations (2.20), respectively. Here, as always, the units are 
such at X = 2tt = 6.28 . . . 

The equations used to compute the curves marked h = are 
foJ = (zVr)- /2 Lr* + 2ir 112 £° r*\tt] , 

J v = 2(*/ir) 1/a f e- it2 dl, 

J±./z 



(2.20) 



'±v7 



where the upper signs are for the shadow side and the lower ones for 
the illuminated side. The computations are made easier by the relations 

(fo^)- - (toJ) + = 2, 
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where the subscripts " — " and "+" denote opposite points on the il- 
luminated side and shadow side, respectively, of the half-plane. These 
relations follow from (2.20).* The radius vector r from the origin is equal 
to — y on the half-plane. Equations (2.20) follow quite readily from (2.3). 
The values of J and J v for an arbitrary angle of incidence and h = 
are given by expressions (6.6) and (6.22), respectively. 

All of the curves in Fig. 2.4, even | faJ \ far h = 1,000, eventually 
approach the value 2 far down on the illuminated side. It may be 
shown that | J v | and | faJ | for the half-plane decrease like {irr)~ m and 
l/(2x 1/2 r 3/2 ), respectively, deep in the shadow. Hence as we go to the 
right in Fig. 2.4, the dashed curves will eventually cross over and lie 
above the solid curves, which decrease exponentially. The larger h, 
the lower and flatter is the curve for | fa J \. 

Now we turn to the diffraction field at points far to the right of the 
cylinder. When | ip | « 1, so that we are not too far from the shadow 
boundary, and h is large, the field is given by (2.10), or by its analogue 
for vertical polarization. In order to get acquainted with (2.10) we first 
examine the field when | i^ [ <3C 1 but \f/ p ^ 1. 

When we are so far behind the cylinder that \f/ 2 p » 1 even though 
I }p | <$C 1, the asymptotic expressions (2.4) show that 

(c(p)/\{/, ^ negative (shadow) 



>S',J 



e~ ix + ^p , yp positive 



(2.21) 



Substitution of (2.21) in (2.10) gives 

c(p)/t 1/3 U(t) + ^ exp (*Y 3 /3), * < 



E ~ < 



cipVi 1 ' 3 U(r) + i) exp (ir 3 /3) + e- ix , + > 



(2.22) 



The presence of c(p) shows that the total wave may be regarded as the 
sum of the incident wave and a wave spreading out from the crest. The 
crest wave is the analogue of the edge wave, c(r)/<p, for the half-plane. 
In fact, when we are in the region where the 1/t in (2.22) is the most 
important term, the crest wave is approximately 

c(p)/f (2.23) 

* They also follow from superposition and consideration of the symmetry of 
the currents produced on the half -plane y < 0, x = when —A is impressed. Here 
A denotes the system of currents which flows in the upper half-plane y> 0, x = 
when the incident wave falls on a complete plane at X = 0. 
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Pig. 2.5 — The amplitude of the crest wave may be obtained from these curves 
and expressions (2.24). Here r = h 1 ! 3 xj, where ^ is small. However p must be large 
enough to make ^ 2 p » 1. 

and this corresponds to a half-plane with its edge at the crest of the 
cylinder, t may be small even though we are considering | \p | to be large 
enough to make (2.21) and (2.22) hold, i.e. large enough to make 
| yp \p l '~ » 1. Indeed, multiplying by h iri gives | r \p m » h 1 " which may 
be achieved for small values of t by making p large enough. 

It follows from (2.22) and its analogue for vertical polarization that 
the amplitudes of the crest waves are 



(hp) 

(vp) 



V-l/Z .1/3 



(2T P r' h"< | *(r) + 1/r |, 
(27rp)- 1/2 h m | *,(r) + 1/r |, 



(2.24) 



where the expression (2.9) for c{p) has been used. The last factors in 

(2.24) may be computed from Table 2.1. They are plotted in Fig. 2.5. 

When we go deep down into the shadow where t is large and negative 
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we see from (2.13) that 

| ¥( T ) 4- 1/r | ~ 2.0:5 
I ¥„(t) + 1/x I — 3.42 . 

so that the absolute value of the field is 



2.025t 



0.882t 



(hp) 
(vp) 



E\ 
H\ 



1/3 



(27rp)" 1/2 ti 13 2.03 exp (2.025 tf'ty), 



i/a 



(27rp)" 1/2 ti' 3 3.42 exp (0.882 tf'ty). 



(2.25) 



(2.26) 



where the angle ^ is negative. Thus, as Artmann has pointed out, the 
field decreases exponentially as we go into the shadow. The larger h is, 
the more rapid is the decrease. This supports the statement made earlier 
that it is darker behind a large cylinder than behind a half -plane. 

Comparison of the expressions for the current density and field 
strength for the shadow regions shows that there is a simple approxi- 
mate relation between them. Near the crest of the cylinder, where x is 
small, the radius of curvature is nearly 2h. Hence the tangent to the 
parabola drawn from the point P (located at (p, \[/) deep in the shadow) 
touches the parabola at T where x is approximately — 2bj/. This is 
shown in Fig. 2.6. Replacing x by —2h\f/ in the expressions (2.19) shows 
that the current density at T is proportional to the field at P as given 
by (2.26). It follows that 



(hp) 
(vp) 



E/UJ 
H/J v \ 



1.41 ti 13 (2irp) 



-1/2 



1.87 /i 1/3 (2ttp) 



\-l/2 



(2.27) 



This leads us to picture the field at P as being produced principally 
by the surface currents around T. The effect of the stronger currents 




IA« 



Fig. 2.6 — The held strength at point P (deep in the shadow) specified by the 
polar coordinates (p, \p) is nearly proportional to the current density at the tangent 
point T specified by the rectangular coordinates (x, y). 
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closer to the crest is perhaps blocked out by the curvature of the cylinder. 
For comparison with the horizontal polarization case Ave note that E 
at (p, \p) for an infinitely long current filament at the origin p = is 
given by 

\E/toI | ~ .5(2xp)- 1/2 (2.28) 

where / is the current carried by the filament and the frequency is 
is such that X = 2ir. There is some difficulty with the picture for vertical 
polarization because the current element at T points directly towards P 
and hence should produce very little field there. This is perhaps as- 
sociated with the fact that the (vp) ratio in (2.27) is smaller than the 
(hp) ratio by approximately the factor /i -1 3 . 

We now leave the shadow region and consider the field at points well 
inside the illuminated region. Fig. 2.5 shows that for large positive 
values of t the amplitude of the crest wave tends to increase with t. 
The asymptotic expressions (2.12) show that when t is large and positive 

| ¥( T ) + 1/t I « I *.(t) + 1/t I ~ (ttt) 1/2 - (^) m h m , (2.29) 

and hence the amplitude of the crest wave deep in the illuminated 
region is, from (2.22) and its analogue, 

(hp) | E - e~ ix | ~ (2wp)~ m (tt#) 1/2 , 

(vp) | H - e~ ix | ~ (2wp)- 112 (xrM) 1 ' 2 . {2 ' Si)) 

Since (2.30) is derived from the general expression (2.10) it is subject 
to the restrictions mentioned just below equation (2.11). In particular 
the angle \f/ should be small (but we must still have ^"p » 1 as assumed 
in (2.22)). When ^ is positive, an application of the laws of geometrical 
optics to determine the reflection from the curved surface of the para- 
bolic cylinder leads to the expressions 



(hp) \E - a"*| ~ [ 
(vp) \H — c~ tx \ ~ 



h tan frfr/2) " 11 ' J 
h tan GA/2)" 11 



sec (iff/2), $ > 

(2.31) 

sec («A/2), ^ > 



for the reflected wave. When \J/ is small these expressions reduce to 
(2.30) as they should. 
Expressions (2.31) may also be obtained from our analysis by start- 

12 In our two-dimensional case the calculation of the required radius of curva- 
ture, etc., is not difficult. General theorems dealing with problems of this sort 
and references to earlier work are given in the paper, A General Divergence 
Formula, H. J. Riblet and C. B. Barker, J. Appl. Phys. 19, pp. 63-70, 1948. 
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ing with equations (7.3G) and (7.56). Furthermore, it may be verified 
that the phase angles of the reflected waves as computed from (7.30) 
and (7.56) agree with those computed from geometrical optics when 
reflection coefficients of —1 and +1 are assumed for hp and vp, re- 
spectively. 

The amplitudes of the crest waves for h = 1,000 and h = are shown 
for hp in Fig. 2.7 and for vp in Fig. 2.8. Of course, when h = the crest 
wave reduces to the edge wave from a half-plane. The curves for h = 
1,000 were computed from equations (2.24) and the curves of Fig. 2.5 
(or their equivalent when t is small). The curves for h = were com- 
puted from 



E | ~ (2tt P )- 



1 + l 



2 sin GA/2) 2 cos (^/2) 



rp < 



(hp) 



e - ,-" I ~ (27r P r i/2 



1 



1 



// 



(2tt P ) 



— 1/2 



2 sin ty/2) 2 cos (tf/2) 



1 



iA > 



(2.32) 



(vp) 



// - e 



(2t P ) 



-1/2 



2 sin W/2) 2 cos (^/2) 



1 



2 sin OA/2) 2 cos ty/2) 



>£ < 



f > 



which follow from (2.1), (2.2), (2.4) and (2.6). 

From equation (2.21) onward we have been discussing the field for 
values of ^ and p such that pip 1 ^> 1. For these values the concept of the 
crest wave is helpful in visualizing the behavior of the field. Now we 
consider the field at points close to the boundary of the geometric shadow 
far behind the cylinder. This is the region in which Artmann 6 was es- 
pecially interested. His results for the shift of the field may be obtained 
from (2.10) and its analogue by taking | \f/ | to be very small. 

At the shadow boundary \p = and [exp (— ix) + S\) P = X A.- Hence 
the region of interest at present is in the neighborhood of the point 
point i\ = 0, | exp ( — ix) + & | = ^-2 °f Fig. 2.2. A magnified view of 
this region showing the shift of the field is given in Fig. 2.9. The figure 
shows that, for a given value of p\f/, \ E \ for hp is less than | H | for vp. 
As Artmann has pointed out, this is to be expected since the reflection 
coefficient for E (hp) is roughly —1 and the reflected wave therefore 
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Fig. 2.9 — Behavior of | E | and I H | on the shadow boundary far behind the 
half-plane or parabolic cylinder. This is a magnified view of the region around 
t { = 0, ! exp (-ix) + S, | = Vi in Fig. 2.2. 

tends to cancel the direct wave when \f/ is very small. On the other hand, 
the reflection coefficient for H (vp) is +1 and the reflected wave tends 
to add to the direct wave. 

The distances 1.71/i 1/3 and -1.49/i 1/3 appearing in Fig. 2.9 are the 
amounts, measured in units for which the wavelength is 2ir, by which 
| E | and | H | are shifted by the curvature of the parabolic cylinder. 
If y' — h' is the shift in meters for | E | and if the radius of curvature 
of the crest is a = 2/i' meters, Fig. 2.9 gives 0(y' - h') = 1.71 W) 1 ' 3 = 
1.71 (/3a/2) 1/3 where = 2tt/X. Thus </ - h' - 0.399X (a/X) 1/3 meters. 
The corresponding shift for | H \ is -0.34GX(a/X) 1/3 meters. Artmann 
gives the values 0.39 and —0.20 for the respective coefficients. The 
discrepancy between —0.346 and —0.20 is cause for worry because it 
seems to indicate either a mistake in our work, which I have been un- 
able to locate, or a shortcoming in the approximations made by Artmann 
for the case of vertical polarization. 

As h approaches zero the parabolic cylinder becomes a half-plane 
and the curves d and c should approach curves b and c, respectively. 
According to Fig. 2.9 both d and e approach curve a. This failure is an 
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indication of the errors introduced by the approximations used in the 
derivation of (2.10) and its analogue. 

3. RADIO PROPOGATION OVER A SUCCESSION OF RIDGES ON THE EARTH'S 

SURFACE 

The results mentioned in Section 1 concerning propagation over a 
succession of ridges may be obtained from the expressions and curves of 
Section 2 as follows: Consider the situation shown in Fig. 3.1. Let a 
radio wave start out from a transmitter at T. We assume that by the 
time it arrives at the first ridge at P it has become equivalent to a 
plane wave of amplitude A/C traveling in the direction TP, where A is 
a constant depending on the strength of the transmitter. For the sake 
of simplicity the waves reflected from the ground are neglected. In a 
more careful study they would have to be included.* 

In order to calculate the strength of the wave at the second ridge, 
we assume it to be a crest wave coming from P. Let G denote the value 
of | E | (we assume the case of horizontal polarization since the reflection 
coefficient of physical materials approaches —1 for almost grazing in- 
cidence) at Q corresponding to a plane wave of unit amplitude incident 
on P. From Fig. 3.1 we see that the values of \J/ and p to be used in com- 
puting G are \J/ « — t/R, p = 2rC/\, X = wavelength, R = radius of 
earth. The value of h depends upon the radius of curvature of the ridge: 
2h = 2t (radius of curvature)/X. 




Fig. 3.1 — Diagram showing ridges at P and Q which diffract the radio wave 
starting from T so that a portion of it is received at S. 



A method for doing this (for one hill) is given in Reference 1, page 417. 
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The amplitude of the wave striking the ridge at Q is AG/(\/2. The 
■\/2 comes from the horizontal sidewise spreading of the wave in going 
from I to 2t. If we were dealing with the energy instead of the amplitude, 
the factor would be 2 instead of \/2. When this wave is assumed to be 
plane and traveling in the PQ direction, similar reasoning shows that 
the amplitude of the disturbance at the receiver S is AG 2 /f\/S. 

If, instead of two ridges at P and Q, as shown in Fig. 3.1, there are N 
ridges between the transmitter at T and the receiver at S, the amplitude 
of the radio wave at S is AG N /C\/N -\- 1. The distance between T and S 
is approximately (N + \)C, and the free space amplitude at S is A/ 
(N + 1)6. Hence 

Actual Amplitude at S nNr-KT i i\i/2 /o n 

— — — -j —^ = (7 (.iv -r- 1; . (.6.1) 

tree Space Amplitude at b 

The actual field at S is therefore 

20 N logio (1/G) - 10 logic (N + 1) (3.2) 

db below the free space field. 

As an example, let us assume a distance of 280 miles between the 
transmitter and receiver, and a distance of 40 miles between successive 
ridges. This gives N = ti. For a wavelength of 10 meters and a radius 
of curvature of 100 meters for the diffracting ridges, the formulas of 
Section 2 show that the ridges behave like half-planes and that G ~ 
0.227. Equation (3.2) then says that, for a distance of 280 miles and a 
wavelength of 10 meters the actual field should be about 69 db below 
the free space field. Although this is in fair agreement with the ex- 
perimental results, calculations for other distances indicate that the 
field strengths predicted by (3.2) tend to be smaller than the ones 
observed. 

When the work is carried through for X = 1 meter and a distance of 
280 miles, (3.2) says that the field is 120 db below free space. The ob- 
served fields are 70 ± 15 db below the free space value. 

These figures suggest that the roughness of the earth's surface might 
possibly account for transmission far beyond the horizon for wavelengths 
of the order of 10 meters. For wavelengths of the order of 1 meter either 
the approximations leading to (3.2) break down or some other explana- 
tion is required. 

4. SERIES FOR THE ELECTROMAGNETIC FIELD 

Here we set down series for the electromagnetic field when a plane 
wave strikes a perfectly conducting parabolic cylinder. Since Epstein's 
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classical work deals with the general case of finite conductivity, the 
series we use are special cases of the ones discussed by him. 

The parabolic coordinates (£, 77) which we shall use are related to the 
rectangular coordinates (x, y) and polar coordinates (r, <p) as follows: 

x + iy - (£ + iiifjli = r e 1 *, 

x = £77 = r cos tp, y = (rf - £ 2 )/2 = r sin V, 
r = (*» + ff* = tf + „ 2 )/2, 
dx 2 + <fy 2 = (f 1 + t, 2 ) (tf£ 2 + dr, 2 ) = 2r{dt + d» 2 ), 

£ = (2r) 1/2 cos^/2 + tt/4), 
77 = (2r) 1/2 sin(*>/2 + tt/4). 

The lines 77 = constant are a series of downward-curving confocal 
parabolas having their focus at the origin. The parabolic cylinder x~ = 
4h(h — y) is given by 77 2 = 2/i. This special value of 77 will be called 770 : 



(4.1) 



770 = C2h) m ^ 0, h = ^/2. 



(4.2) 



When 770 = 0, the cylinder reduces to the half-plane x = 0, y ^ 0. The 
lines £ = constant are halves of upward-curving confocal parabolas 
having their common focus at the origin. Outside the cylinder 77 > 770 ^ 0, 
so 77 is always positive in our work. £ is positive in the half-plane x > 
and negative in x < 0. 

For much of our work we shall assume the incident wave to come in 
at the angle 0, ^ 6 ^t, as shown in Fig. 4.1. As mentioned in Section 



4= CONSTANT 



/7 = CONSTANT 




V = Vo 



Fig. 4.1 — This diagram shows the angle of the incident wave and the surface 
of the perfectly conducting C3'linder x 2 = ih(h — y) (or 77 = r;o). 
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2, the field quantities are assumed to depend upon the time through the 
factor exp(i'ut) where co is the radian frequency. 

The wave equations for horizontal and vertical polarization are, 
respectively, 

S + 8 + ( * 2 + " )E = ° a3) 

^ + $? + (? 2 + 1*) A = (4-4) 

where, as explained in Section 2, the unit of length has been chosen so 
that the wavelength X = 2t. On the surface of the perfectly conducting 
cylinder, i.e. for 77 = rj , we must have E = and dH/dij = 0. When 
E and H are known the remaining components of the field may be 
computed from Maxwell's equations. 
Special solutions of (4.3) (and (4.4)) are 

exp [itf - if)/2] U n fa w ) tf»(*r 1/2 ), (4.5) 

exp [»(*' - ?)/2] U n (& m ) W n ( v i~ m ), (4.6) 

where i m stands for exp (tV/4) and £/„(z), T7 n (2) satisfy the equation 

^W _ 22 ^M + 2nr„(«) = 0. (4.7) 

dz- dz 

Another solution of (4.7) with which we shall be concerned is V n (z). 
These three solutions are defined by contour integrals of the form 

(27rr)~ 1 f exp [f(t)) dt where /(/) = -t 2 + 2zt - (n + 1) log «. 

The path of integration for t/ n (z) comes in from - *> where arg £ = —v, 
encircles the origin counterclockwise and runs out to — <*> with arg 
t = ir. The path for V„(z) runs from — °o where arg t = x to + 00 where 
arg £ = 0, and the path for W n (z) runs from + » to - *> where arg 
i = — 7T. The integrals are written at greater length in equations (9.1) 
and the paths of integration are shown in Fig. 9.1. Since the paths may 
be joined to form a closed path containing no singularities of the inte- 
grand it follows that 

U n (z) + V n (z) + W n (z) = 0. (4.8) 

When n is a non-negative integer 

UM = H n (z)/n! = ^" e' 2 f- e*' (4.9) 

n! dz n 
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where H n {z) is Hermite's polynomial. When z becomes very large the 
leading terms in (9.17) and (9.16) give 

U n {z) ~ (2z) n /n>, (4.10) 

W n ( V i- 112 ) ~ itf'VvT* 1 e-^Ar 1 ' 2 . (4.11) 

In order to obtain a series for the incident wave 
exp [— ix sin 8 + iy cos 6] 
shown in Fig. 4.1 we consider the special case = 0. In this case the wave 
is simply exp(iy) or exp [i(-q — £ 2 )/2] and may be obtained by setting 
n = in (4.5). This suggests that the incident wave may be expressed 
as the sum of terms like (4.5). The series turns out to be 

exp [ — ix sin + iy cos 0] 

= exp [ — ££ij sin + i cos 0(ij 2 — £ 2 )/2] 

= exp [- izz' sin - cos 6(z 2 + z ,2 )/2] (4.12) 



where 



= e^sec (6/2) X n!(- iw/2) n U n {z)U n (z') 

w = tan {6/2), 

z =# 1/2 , (4.13) 

z' = nf 1 *. 

This series has been studied by a number of writers. It goes back to 
Mehler 13 who obtained it by evaluating the integral 



tT 1 ' 2 e xi ( M exp [- (t - iy) 2 - (x + iat) 2 ] 

J— 00 



dt 



first in closed form, and then as a series (by using the generating func- 
tion exp [— {—iat) 2 + 2( — iat)x] for H n (x) and integrating term wise). 
This leads to 



(1 - aT 112 exp 



2xya — (x 2 + y 2 )a 
1 - a 2 



(4-14) 
= £ H n (x)H n (y) a n /2 n n! 



which is equivalent to (4.12). Since (4.14) converges when \a\ < 1, 
(4.12) converges when | wj | < 1 or | 6 \ < ir/2. 

13 Reihenentwicklungen nach Laplaceschen Functionen hoher Ordnung, J. 
Reine Angew., Math., 66, pp. 161-176, 1866. 
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When the incident wave strikes the cylinder the reflected wave has 
some of the characteristics of a wave spreading radially outwards. Such 
a wave contains the factor exp ( - ir) = exp [ - i(g + t? 2 )/2] . Consideration 
of the exponential factors in (4.6) and (4.11) suggests that the reflected 
wave may he expressed as the sum of terms of the form (4.G). The co- 
efficients in this series are to be determined so that E = or dH/di) = 
at the surface r\ = t? , the incident wave being represented by (4.12). 

For the case of horizontal polarization this procedure gives 

00 

E = e iv sec (6/2) £ n!(- iw/2) n U n (z) [U n (z') (4 15) 

-W,Xz')U n (z' )/W n (z' )] 
E = exp [ — ix sin + iy cos 0] 

- e tv sec (0/2) Y, n/(- iw/2) n U n {z)W n {z')U n (z,)/W n {z ) (4 ' 16) 



for the complete field. These are special cases of Epstein's results. Here 
z'o is the value of z' which corresponds to the surface of the cylinder: 

z ' = r% = (2h/i) m (4.17) 

The entries for regions II and IF (these are regions in the ra-plane 
(m = n + 1) which, as Figs. 11.2 and 12.2 show, contain the large posi- 
tive values of n) in Tables 12.2 and 12.4 of Section 12 may be used to 
show that as n — > <» 

V n (z' )/W n (z' ) ~ r 2n exp [2 7?0 (2m) 1/2 ], 

U n (z' )/W n U) = -1 - V n (z' )/W n (z' ), 

UMWM - - exp [ " iy ~ v(2n/i)1,i] 



(4.18) 



Since 



4[I\1 + n/2)]* 
{exp[- £(2/iA) 1/2 ] + i- n exp [$(2n/i) 112 ]}. 

n//fr(l + n/2)] 2 ~ 2" {im/2T m 



the series in (4.15) and (4.16) converge if | w \ = \ tan 0/2 \ < 1. 

Series for H similar to those of (4.15) and (4.16) may be obtained for 
the case of vertical polarization. The boundary condition at t] = ij 
is now dH/d-q = 0. It is convenient to introduce the functions 'U„(z), 
'V n (z), 'W n (z) defined by 

f U n (z) = -zU n (z) + dU n (z)/dz (4.19) 
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and the two other equations obtained when U is replaced by V and W. 
The prime is placed in front of U instead of behind to avoid mistaking 
'U„(z) for dll n (z)/dz. The function 'U n (z') makes its appearance when 
dH/dri is calculated for the boundary condition. Since iy = — (z 2 + z /2 )/2 
we have 

^ e iv U n (z') = C m e iv 'U n (z'). (4.20) 

077 

The analogues of (4.15) and (4.16) for vertical polarization are (as- 
suming now that H for the incident wave is of unit amplitude). 

// = c iu sec (0/2) £ n!(-iw/2) n U n (z) (4 21) 

[U n (z') - W n (z')'U n (zo)/'W n (z' )}, 

= exp [—ix sin 6 + iy cos 8] 

(4.22) 
- e''"sec {d/2)Y,n!{-iwm n U n {z)W n {z')'U n {z,)/'W n {z a ), 


and these series converge if \w\ < 1. 

If the parabolic cylinder is merely a good conductor, instead of being 
perfect, the boundary conditions at 77 = 770 are approximately E = 
— f#f, E{. = £#. Here E^ and H| denote the £ components of the elec- 
tric and magnetic intensities and f is the intrinsic impedance 

f = [*W(0 + *««)] 1/2 (4.23) 

of the cylinder material, f is assumed to be small compared to the 
intrinsic impedance U = (moAo)" 2 = wmo (since X = 2tt) of the external 
medium. In these expressions n, e, g are the permeability, dielectric 
constant, and conductivity of the cylinder; and (i and e refer to the 
external medium. 
When we set 

(T = l -^4-r 7 o 2 ) 1/2 fo/f, (4.24) 

r = t- I/2 (r ! + .o 2 ) 1/2 r/ro, 

the boundary conditon for hp becomes <jE = —dE/dz' at z' = z' Q . When 
a is assumed to be constant we obtain 

E = a* sec(0/2) fj n/(-tto/2)"l7.« 

(4.25) 

|M - TF„(sO[<r^«(z5) + 'UJ&WWJjA + 'F.Cfc)]! 



" Electromagnetic Waves, S. A. Schelkunoff, D. Von Nostrand Co., N. Y. (1943) 
p. 89. See also G. A. Hufford, Quart. Appl. Math. 9, pp. 391-403, 1952, where ref- 
erence is made to the work of Leontovich and Fock. 
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which reduces to (4.15) when a — > go . The constancy of <x may be achieved 
by either taking the properties of the cylinder material to change in a 
suitable way or, roughly, by taking 770 (and hence h) to be so large that 
only the nearly constant values of a at the crest of the cylinder have an 
effect on the result (assuming 9 = 7r/2, and restricting our attention 
to the region near the shadow boundary). 

The corresponding expression for vertical polarization may be ob- 
tained from (4.25) by replacing E and <r by // and t, respectively. We 
shall refer to (4.25) and its analogue later in connection with the field 
in the shadow (Section 7) and with Fock's investigation of the surface 
currents on gently curved conductors (Section 6). 

5. INTEGRALS FOR THE FIELD 

When the curvature of the cylinder is small, i.e., when h is many 
wavelengths, the series of Section 4 converge slowly. The work initiated 
by G. N. Watson 1 on the smooth earth problem suggests that we con- 
vert the series into contour integrals with n as the complex variable of 
integration. When this is done we get an integral with the path of in- 
tegration L\ shown in Fig. 5.1. Thus, for example, expression (4.16) for 
E is transformed into 



E = exp [—ix sin -f- iy cos 0] 





sec- 



2% 



[ ( l 4) n ^^^ U n (z)W n (zW n (zo) dn/W n (z' ). 
jli \ 2 / sin 7m 



(5.1) 



At first sight it seems that not much can be done with this integral 
because the integral obtained by deforming Li into hi does not converge 
(this is explained in the discussion of Table 5.1). However, some ex- 



L 3 


L 2 


( 




n-PLANE 

L, 








-2 -J) 


} 


2 « 
n = -lh = 


2 

-*«,72 






- 


r^ZEROS OF W n (Zq) 



Fig. 5.1 — Paths of integration in the complex n-plane. 



16 Proc. Roy. Soc, London (A) 95, p. 83, 1918. 
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perimentation shows that if we set' U n {z ) = -V n (z' ) -W„(z' ) in (4.16), 
the series splits into two series, one of which may be summed and the 
other may be converted into an integral along the path L 2 of Fig. 5.1: 

E = exp [-ix sin + iy cos 0] + & + 82(h), 

Si = e iu sec (0/2) Z n!(-iw/2)"U n (z)W n (z'), (52) 

So_(h) = c/'sec (0/2) Y J n!{-iw/2) n lJ n {z)W n {z')V n {z )/W n {.z' (i ). 

The series for & may be summed by replacing U n (z) and W n (z') by 
their expressions (9.19) in terms of definite integrals, and interchanging 
the order of summation and integration. The resulting series may be 
summed and the integrations performed. The result is Sommerfeld's 
integral for a diffracted plane wave: 



8l = _ (l - A) x/v* -'+*«>■' r e- ita dt, 

T 1 =,cos^-^in-^=(2.)-sin(^4-I 



(5.3) 



The inversion of the order of summation and integration may be 
justified when | w \ = \ tan(0/2) | < 1 (in which case the series in (5.2) 
converge) by using (1) the result that | R N \ < \ a N /N! | when a is real in 

n! 
and (2) the inequality 

ft N ex V [-t a + 2 ll2 bt]diT:- 

< f e -^-**r dt f exp [-2«/ 2 + 2 m bt] dt 

Jo J— <*> 

< A2- N b- m N- vi N! exp(2bN 112 ) 

This inequality holds when N » b 2 and at the same time N » 1. The 
value of A is independent of N and b is a number which exceeds | £ |. 
In this work the parameter a has been arbitrarily introduced; and has 
then been chosen so as to make the product of the two integrals a mini- 
mum when iV is large. This value of a is b N ". 

When the series (5.2) for S 2 (h) is converted into a contour integral 
taken along the path L u by the procedure used to obtain (5.1), it is 
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seen that L\ may be deformed into L-> and we obtain 



sm . » /" m"Efe±«i7. W ir.(rtF.&0*^r.(4) (5 ' 4) ' 

2z Jt, \ 2 / sin 7m 

Whether a particular integrand, such as the one shown in (5.4), con- 
verges at the ends n = ±i°° of L 2 can often be decided from Table 5.1. 
This table gives a rough idea of the behavior of the various functions in 
terms of powers of i. For example, if the integrand should turn out to 
be proportional to i" = exp(iirn/2) at n = i<x>, the integral will con- 
verge like exp ( — t \ n | /2). 



Table 5.1 





Order of Magnitude — Rough Approximation 


Function 






near » = «'« 


near » = — i» 


i n 





co 


i-» 


CO 





sin irn 


i~ 2n 


j2n 


r(n + 1) 


i n 


i~" 


U n (z) 


j-3»/S 


J3nl2 


V n (z) 


7 "-3n/2 


i-nl2 


W„(z) 


j-n/2 


V3n/2 



The approximations for T(n -\- 1) follow from its asymptotic expres- 
sion, and those for the parabolic cylinder functions come from Tables 
12.2 and 12.4. The entries for the cylinder functions may also be sur- 
mised from expressions (9.4) which hold for z = 0. 

Table 5.1 may be used to show that the integrand in (5.1) is of the 
order of i" as n — > — i « . Hence there is no hope of deforming Li into 
L 2 in this case. On the other hand, the integrand in (5.4) is of the order 



of i" as n —> ix> and of i " as n 



— ICC 



and therefore (5.4) converges 



exponentially. In fact, it converges for all real positive values of w = 
tan 0/2. This enables us to obtain an expression for the field which 
holds for < 6 < t (i.e. it is not subject to the restriction | w \ < 1 re- 
quired by (4.16)). This expression, which is fundamental for our work, 
has the form 



E = exp [—ix sin 6 + iy cos 6] + Si + S- 2 (h). 



(5.5) 



Here Si and S?(h) are given by (5.3) and (5.4), respectively. 

In working with (5.5) it is sometimes convenient to use the expression 
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exp [—ix sin 4- iy cos 0] + Si 

f Ti _•„ ( 5 -6) 
= (r*/ir) 1/2 exp [ — m; sin + iy cos 0] / <? ' dt 

J— 00 

which follows from (5.3) and 

'expC-iflA = Gr/0 1 ". (5.7) 



£ 



The development leading to (5.5) shows that is satisfies the boundary 
condition E = at j? = t/ f or < w < 1. That (5.5) also satisfies the 
condition for the extended range < w < «> follows immediately from 



e sec - . / . \ n «/ i t, 
2 f /iw\ T(n -f- 1) 



J L ,\2/ sii 



U n (z)V n (J) dn 



2i Jl 2 \2 / sin 7m 

- ri o (5-8) 

= -(i/ir) 11 ' 2 exp [ — iz sin + % cos 0] / exp (—ft*) (ft 

J— 00 

= —exp [— ix sin + iy cos 0] — Si 

when we note that setting z' = z'o reduces S 2 (h) to the left hand side of 
(5.8) (with z' = z'o). 

Equation (5.8) is due to T. M. Cherry 16 who obtained it by expressing 
the cylinder functions as integrals and interchanging the order of in- 
gration (he works with the function D n (z) of our equations (9.2)). Sub- 
stituting the integrals (9.19) for U„(z) and V n (z') in (5.8) and inter- 
changing the order of integration leads to a similar derivation. Equation 
(5.8) may also be obtained by deforming L 2 into Li when < w < 1 
and into L 3 when 1 < w < <x> . This leads to the two series 

e iv sec 6 - 2 (-iw/2) n n!U„(z)V n (z'), (5.9) 

2 »i=o 

-e*sec| Z (-iw/2) n [T(n + l)U n (z)]V n (z') (5.10) 

2 n=-l 

which may be summed in much the same way as was (5.2) for Si. 

An expression for E which is useful in the study of the current density 
on the surface of the cylinder may be obtained from (5.5) by combining 
expression (5.8) for exp [—ix sin + iy cos 0] 4- Si with expression (5.4) 

18 Expansions in Terms of Parabolic Cylinder Functions, Proc. Edinburgh 
Math. Soc, Ser. 2, 8, pp. 50-65, 1948. 



DIFFRACTION OF RADIO WAVES BY A PARABOLIC CYLINDER 449 

for 82(h): 



JL(t 



2i Jl., \ 2 / sin 7m 



r y w (gp) v,(g ; ) 1 ,, 

lW n (z'o) W n &)\ 



(5.11) 



When w exceeds unity (or when w = 1 and £ > ?? ^ j» ^ 0) in (5.11), 

it may be verified with the help of Tables 12.2 and 12.4 that Li may be 
deformed into L 3 + L^. When n is a negative integer the quantity within 
the brackets in (5.11) vanishes because of (4.8) and because U n (z) = 0. 
The contribution of L 3 is zero since it encloses no poles. The contribution 
of La is equal to the sum of the residues at the poles given by W„(zo) = 0. 
Hence, when w > 1, 



E = —jre iy sec - 2Z 

2 8 = 1 



iwY r(n + l)U n (z)W n (z')V n (z' )' 



2 / sin 7mdPT n (2o)/dn 



(5.12) 



where n = n, is the sth zero of W n (z ). This series also converges when 
w — 1 and £ > rj ^ 770 ^ (which is roughly the shadow region). The 
preceding inequality does not necessarily specify the complete region 
of convergence. 

Cherry 16 has also pointed out that the expression for a plane wave 
given by A. Erdelyi 1 , namely (in our notation) 

exp [ — ix sin 6 + iy cos 9] 

= -llp( (%yrJ*±n luM r.w (513) 

2% Jlo \ 2 / sin irn 

+ U n (-z)V n (-z')]dn, 

may be regarded as the sum of the negative of (5.8) and a similar ex 
pression with £ and 77 replaced by — £ and —77. In informal discussions 
with the writer, Prof. Erdelyi has pointed out that the work leading to 
our expression (5.5) for the field may be considerably shortened by 
starting with some known integral for the impressed field, such as (5.13) 
or a related result. One way of doing this is to take 

exp [—ix sin + iy cos 6] + Si, 



> 7 Proc. Roy. Soc Edinburgh, 61, pp. 61-70, 1941. 
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as given by the left hand side of (5.8), to be the impressed field in the 
equation 

E = impressed field 4- reflected field 

From the form of (5.8) and the discussion of expression (4.6) (given 
between equations (4.14) and (4.15)) we are led to assume the reflected 
field to be an outgoing wave of the form 



e sec - 
2i 



f (^\ ^±^ U n (z)W n (z')a(n) dn 
Jl 2 \ 2 / sin irn 



where a(n) must be chosen so as to make E vanish on the surface of the 
cylinder. This gives a(n) = V n (za)/W„(z'o) and leads directly to the 
expression (5.11) for E. 

When the incident wave is vertically polarized, integrals for H may 
be obtained from the series of Section 4 in much the same manner as 
were the integrals for E. The analogues of the earlier results are 

H = exp(— ix sin + iij cos 0) 



e rv sec - 



I /toV r£+0 VJlAWM yvJA dn/'WM), 

J/,! \ 2 / sin irn 



(5.14) 



% 

II = exp (-ix sin + iy cos 0) + Si + S 3 (h), (5.15) 



'7„(4) dn/'W n (z'o), 



e" sec - 



2t Ji,,\2/ sm Tfi (5.17) 

[ ^n^) _ Vn(z') 1 
i'Wniz'o) W n { Z ')\ dn > 



- r(tU»/2)"r(n + l)Unb)WnWVn(&' 



18) 



# = -Tre'" sec (0/2) £ : _, w , n ,, (5. 

In these formulas 'U n (z), etc. are defined by (4.19); w, z, z' by (4.13), 
z'o by (4.17). In (5.14) w is restricted to < w < 1. In (5.15) 5 8 (A) is 
given by (5.16), and w may be anywhere in < w < °o. In (5.17) w 
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may also lie anywhere in < w < », but in (5.18) it is restricted to 
1 < w < oo except when £ > rj (roughly the shadow region) in which 
case w may be unity. In (5.18) n = n' s is the sth zero of 'JF n (zo)." The 
zeros of 'W„(z a ) interlace those of W n (z ) shown in Fig. 5.1. 

When /; = ijo/2 = the parabolic cylinder degenerates into a half- 
plane and our solutions reduce to Sommerfeld's expressions for waves 
diffracted by a half-plane. It may be shown that if 



r,-,o»| + €m|-W w «ni(!^+i+») I 



(5.19) 



we have 



a(0) = -{i/w) 1 '* exp [ix sin 6 4- iy cos d] [ c~" 2 dt, 

Jt * (5.20) 

&(0) = -S,(0). 

When this expression for So(0) is added to (5.6) Ave obtain Sommerfeld's 
result for the case of horizontal polarization. 

One may verify that the series (5.12) leads to Sommerfeld's result as 
Zo approaches zero. By neglecting 0(z 2 ) terms in (9.3) and setting n B = 
— 2s + V» f° r the sth zero of W n (z ) we may obtain the following rela- 
tions which are needed in the course of the verification ' - 

p. = -4&oT(« 4- l/2)Arr($) + • ■ • , 

dW„(z' )/dn at n. = T(s)/4 + • • • , 

Yn(z' a ) at n 8 = -2iz' T(s 4- 1/2)/*- 4- • • • , (5.21) 

V n (z' ) 



r vm I 

[sin irndW„(z' )/dn} n= „ s 



= 2/tt4- 



6. SURFACE CURRENTS ON THE CYLINDER 

As shown in Fig. 6.1, the surface current J on the perfectly conducting 
cylinder 77 = 770 is parallel to the crest of the cylinder (and to the elec- 
tric intensity E) when the incident wave is horizontally polarized. We 
have from Maxwell's equations in parabolic coordinates 

J = [-HJn, = (tfo)" 1 (2r)- 1/2 [dE/dv]^. (6.1) 

Here H^ is the component of magnetic intensity in the f -direction, 
fn is the intrinsic impedance of free space given by fo = (moAo) 1/2 ™ WjUo 
where the second part of the equation follows from 2ir/X = u(mo6o) x/s 
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and X = 2tt. The E and the H of the incident wave are related by 
H Q = Eo/to = 1/fo since E = 1. In rational MKS units f = 120tt ohms. 
The derivative in (6.1) may be obtained by differentiating expression 
(5.11) for E and then setting ?? = vo- Use of the Wronskian (9.9) for 
V n (z' ), W n {z' ) then takes (6.1) into 

r „j = M [ -^— (iw) n U n (z)/W n (z' ) (6.2) 

Jl.j sin irn 

where w = tan (6/2), L 2 is showm in Fig. 5.1, 2 and z' a are given by (4.13) 
and (4.17), and 

Mo = (i/Sirr) m e- ir sec | , 

r = (r + 7, 2 )/2. 

In this Section r will be restricted to mean a radius vector drawn to the 
trace of the cylinder on the (re, y) plane of Fig. 6.1. 

Closing Li on the right and on the left gives the two series 



£ J = 2iM E (-iw) n U n (z)/W n (z' ), < w < 1, 

n=0 

(iw)"U n (z) 



f J = -2iwM £ 



«=1 - 



sin tt?i 3W„(2o)/dn 



] . 



(6.4) 
1 < to < oo, (6.5) 



where n 8 is the sth zero of W n (z' ) regarded as a function of n. 
For the half-plane case j? = 0, and (9.4) gives 

W n (0) = -i n /2T(l + nil). 

In this case the series (6.4) may be expressed as an infinite integral 




Fig. 6.1 — Relationship between surface current density J and electromagnetic 
field when incident wave is horizontally polarized. E and J are normal to the 
plane of the paper. 



DIFFRACTION OF RADIO WAVES BY A PARABOLIC CYLINDER 453 

when the integral for r(l -f w/2) is inserted and the sum (9.22) [i.e., 
the sum for the generating function of U„(z)] used. Integrating part of 
the result gives 

fo-7 = (2/hr) i;2 cos^[e- <r 

' ' (6.6) 



- 2# sin i «"* cose T exp ( - U*) dt\ . 

2 JZ sin (6/1) J 



In (6.6) r is the distance along the half-plane as measured from the 
edge: r = £ 2 /2 = \ y \. Positive values of £ correspond to the shadowed 
side of the half-plane and negative values to the illuminated side. With 
this interpretation (6.6) agrees with the current density obtained from 
Sommerfeld's expression for the field. 

Although (6.6) has been derived from (6.2) on the assumption that 
< w < 1 it also holds f or < w < « as may be shown by analytic 
continuation. Again, (6.6) may be obtained from (6.5). 

Since the factor r~ 1/2 comes from the multiplier M in (6.4), it is pos- 
sible that (6.6) may give one an idea of how the current density behaves 
near the crest of a thin cylinder which is almost, but not quite, a half- 
plane. Of course, r would have to be interpreted as shown in Fig. 6.1. 

In order to study J when the radius of curvature of the cylinder is 
large compared to a wavelength we consider the case 8 = 7r/2, i.e. w = 1, 
in which the incident wave comes in horizontally. In this case most of 
the variation of the current density occurs near the crest of the cylinder 
where, as it turns out, £ is of the order of t? 1/3 , tjo being large. 

At the beginning of the investigation rough calculations of the inte- 
grand in (6.2), based on the asymptotic expressions of Section 12, sug- 
gested that for small £ and large 7j : 

(a) Most of the contribution to the integral (6.2) comes from the 
neighborhood around point C shown on Fig. 6.2 where m = n + 1 = 
zq 2 /2 = —itfa/2 = —ih. Point C is a critical point associated with the 
asymptotic behavior of W„(z' ). 

(b) The path of steepest descent for (6.2) roughly corresponds to the 
line ACD of Fig. 6.2, C being the high point of the path. Along this 
path Im [f(/o) - f(ti)] = where /(0 = -t 2 4- 2zt - m log t,m = 7i + 1, 
is the function entering the definition (10.1) of the parabolic cylinder 
functions, and f , k are the saddle points of exp [/(/)] . This path in the 
n-plane separates the regions in which \V n (z' ) has different asymptotic 
forms. It is the boundary of region III' in Fig. 12.2 and has been studied 
in Sections 11 and 12. 

Once (6) is verified the truth of (a) follows almost immediately since 
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the path of integration L 2 may be deformed into ACD without passing 
over any singularities of the integrand of (6.2). 

In order to verify (b), we note that the entries in Table 12.3 for W n (jio) 
show that along ACD 

W n (4)~A' ttexp[f(U>)). (6.7) 

Here the expressions for W„(z' ) along ACD are taken to be those cor- 
responding to the regions / b and II shown in Fig. 12.2. In making the 
last approximation in (6.7) we have neglected the slowly varying co- 
efficient of the exponential function in the expression (12.9) for A . 
Since | £ | <5C r} we set £ = in U„(z). Then upon using the values (9.4) 
for U n (0), (6.7) for W n (z' ), and unity for w, we see that the integrand 
of (6.2) behaves like 

i n exp [-/(/ )] ( 68 n 



2 sin (7ro/2)r(l + n/2) ' 

On ACD we have, in dealing with W„(z'„), —St/2 < arg m ^ — t/2. 
Hence, from h + t\ = i~ Ui r) and from f(ta) + f(ti) as calculated from 
(12.9), we have 

exp [/(/„)] - exp (*[/«o) + f(Q] + M/(4) - f(ti)}) 

~ 2" +1 (27r)- 1/2 r(-n/2) exp (- 2? 4- §[/(/„) - /fo)]> 

where we have used the second of expressions (12.10) to evaluate 
exp [m(l - log (m/2))/21. Substitution of (6.9) in (6.8) shows that the 
integrand behaves roughly like 

exp (^ - |[/(f ) - /(/i)]) • (6.10) 

The truth of statement (b) then follows from the fact that the lines o' 
steepest descent of (6.10) in the n-plane are given by Im [f(t ) — f(ti)] = 
0. To see that C is the high point of ACD we use (12.9) to show that 
near C we have 

/(/„) - f(h) W (2/3zo) (zo 2 - 2m) 3 ' 2 

where m = n -\- I. Consequently, /(f ) — j(t\) is real and positive on 
AC [where, near C, arg (z 2 — 2m) = — tt/6] and on CD [where arg 
(z —2m) = —Sir/2, m being in region II' according to the convention 
used in (6.7)]. That C is the high point now follows from (6.10). 

In accordance with statement (a), we must study the form assumed by 
the integrand of (6.2) when n is near point C. 
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When (1) n is near C, (2) 770 is large, and (3) | £ | « 170 we have for the 
various terms in (6.2) 



/"/sin irn ~ 
r(l 4- n/2)W n (z'o) - 

2r(l + n/2)U„(z) ~ i 



■2i l 



U/4) i;3 (2^) 1/ V' 6 exp [-ivo/2}Ai(a), 



exp b 



2 *\i / 6 \2m 



(6.11) 



(6.12) 



where Ai (a) denotes the Airy integral denned by (13.12), m = n + 1, 
arg m is near —ir/2, and 

a = (2/i v 2 ) m (m + ivl/2), o (6.13) 

da = {2/n\) m dn. 

Expression (6.11) comes from (13.21) and expression (6.12) comes 
from region la of Table 12.2 (strictly speaking, region lb should be 
used but point C is so close to arg m = -ir/2 that the simpler expression 
for la may be used). In obtaining (6.12) it is necessary to use the terms 
shown in the expansions (12.5) of / and log t x /U. It may be shown that 
(6.12) also holds for negative values of £. 

We now set w = 1 in (6.2) and change the variable of integration 
from n to a. Substituting for m in (6.12) its expression in terms of a, 
expanding in powers of a and neglecting higher order terms, converts 
the argument of the exponential function into 



if/2 — ^770 — *7 3 /3 + ayi 



1/3 



(6.14) 



B L0 ° 




C is at n = -1-lh 



D -Loo 
Fig. 6.2 — Paths of integration used in studying the current density and diffrac- 
tion pattern when h is large. Path BCD is equivalent to path L 2 of Fig. 5.1. AC 
and CD are boundary lines which mark a change in the asymptotic behavior of 
Wn(zlt). Far out towards .4 the line AC tends to become parallel to BC. 
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where 

7 = £/(2t7o) 1/3 = x/2h w (6.15) 

When our approximations are set in (6.2) we obtain 
UJ » (i/a*)"**- 1 exp [-ifo - iy 3 /3] 

-ooexp (-»2ir/3) (6.16) 

/ exp (i 1 3 ya) da/Ai(a). 

Joo exp (i2ir/3) 

Here we have taken the path of integration in the complex a-plane to be 
the transformed version of the path of steepest descent in the n-plane. 
That the path for (6.16) is still the one of steepest descent when 7 = 
and a is large follows from the asymptotic expansion (13.19) for Ai{a). 
It is interesting to note that near the crest of the cylinder £770 + y 3 /S is 
approximately the distance along the cylinder as measured from the 
crest. 

The expression (2.14) for fo«/ is obtained from (6.16) when Ai{a) is 
transformed by the relations (13.17). 

When 7 is positive the path of integration for a in (6.16) may be closed 
on the left to obtain a convergent series, the terms of which arise from 
the zeros of Ai(a). These zeros lie on the negative real a-axis starting 
with a - ai = —2.338 • • • , a 2 = —4.087 • ■ • The first fifty values of a„ 
and the corresponding values of the derivative Ai (a 8 ) have been tabu- 
lated*. Thus, when 7 is positive, 

°° / -1/3 \ 

WJ « (2/^W'^-^ 3 E eXP /* 7 r (6.17) 

,=1 At (a.) 

The leading term in this series leads to the approximation (2.19) when 
we use Ai f (fli) = .701 • • • Expression (6.17) gives the form assumed by 
'6.5) when w = 1 and h —» °o . For large values of s* 

a, [3tt(4s- l)/8] 2/3 

(6.18) 
Ai'ifl.) (-)V- 1/2 (-a s ) 1/4 

When 7 is large and negative an asymptotic expansion for £ J may 
be obtained by setting the asymptotic expansion (13.19) in (6.16) and 
using the method of steepest descent. It is found that the saddle point 
is at a = a = yi 3 and the slope of the path of the steepest decent 
through it is given by arg (da) = — 57r/12. This leads to the expression 
for ft-/ which appears in (2.17). 

When the incident wave is vertically polarized the magnetic intensity 

* Reference 11, page 424. 
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H is parallel to the crest of the parabolic cylinder. Since the cylinder 
is a perfect conductor, the current density J v on the surface 77 = 170 
is equal in magnitude to H and its direction is that of increasing £. 
Thus setting z' = z' in expression (5.17) for H and using the Wronskian 
(9.9) gives 



J v = N f -J^- (iw) n U n (z)/'W n (zo), 
Jt„ sin Tti 

N = e~ ir (sec °-\ / 2t v \ r = fei 4- £ 2 )/2, 

where 'Fn(«o) is defined by (4.19). 

Closing L 2 on the right and left leads to the analogues of (6.4) and 
(6.5): 

CO 

./„ = 2/JVX (-iw) n U n (z)/'W(zo). < 10 < 1, (6.20) 

71=0 

(t») m W.« 



J„ = -2i7rA r £ 



1 < w < oo, (6.21) 



sin imd'W n (zo)/dn_ 

where n'„ is the sth zero of 'W n (z' ). The zeros of both 'W n (z' ) and PT B (2o) 
are enclosed by the path of integration L 4 shown in Fig. 5.1. 

The current density on a half-plane is obtained by setting z' = in 
(6.20), using 'TP n (0) = W n '(0) = -i"' 1 /Yin/2 4- y 2 ), from (9.4), 
and the generating function series (9.22) for U n (z): 

J v = 2{i/irf n e- iTCOe6 f e~ ili dt (6.22) 

J l sin (9/2) 

where r has the same significance as in (6.6). This agrees with the ex- 
pression obtained from Sommerf eld's result for the half-plane. 

When w = 1 and h is large, the path of steepest descent for (6.19) 
becomes the same as that for the case of horizontal polarization, namely 
ACD of Fig. 6.2. This follows from the fact that, as may be seen from 
(12.2), the controlling exponential functions for 'W n (z' ) and W„(z' Q ) 
are the same. The analogue of (6.16) is obtained by using the approxi- 
mation (13.24) for 'W n (z' ): 

J v pa (l/2iri) exp [— i&ia - iy 3 /3] 

^ exp (-i2W3) ( 6 .23) 

exp (i m ya) da/Ai'(a) 

Jco exp (j2r/3) 

where Ai'ioi) denotes dAi(a)/da and 7 is given by (6.15). 
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For positive values of y (6.23) and d"Ai(a)/d a = aAi(a) lead to 

oo t "1/3 '\ 

J v « -exp l-?^ - Z7 /3J 2- '../ '\ t6.24) 

where a,' = -1.019, a% = -3.248, • ••, etc. are the zeros of Ai'(a).* 
When we use Ai{a[) = 0.5357 the leading term in (6.24) gives (2.19). 
For large values of s* 

a, [3tt(4s - 3)/8] 2/3 , 

Ai{a' s )~ -(-)\- a ' s )- m TT m . 

The expression (2.17) for J v when y is large and positive may be ob- 
tained by applying the method of steepest descent to (6.23). The asymp- 
totic expression for Ai'(a) is obtained by differentiating (13.19). 

When the cylinder is a good conductor, but not perfect, the expression 
for H analogous to (4.25) leads to an integral for J v , obtained from (6.23) 
by substituting Ai'(a) + (Ai(a) for Ai'(a), which is equivalent to one 
given by Fock.f Here f — — (ih) 1 $/{"o is assumed to be small compared 
to unity and f /£ o is the ratio of the intrinsic impedance of the cylin- 
der material to that of free space. Horizontal incidence, = ir/2, is as- 
sumed. 

The analogue of (4.25) for H has the same form as (4.21) except that 
now 'U„(z' ) is replaced by 'U„(z' ) + tU„(z' ), etc. The development 
leading from (4.21) through (5.15), (5.17), (6.19) to (6.23) may be car- 
ried out just as before. The work is also related to that given at the end of 
Section 7 where the effect of finite conductivity on the diffracted wave 
is briefly discussed. 

A series corresponding to (6.24) may be derived from the integral. 
The exponential terms in this series are approximately exp [i m y(a' 8 
— ■C/a'g)], and are similar to those in (7.63). Since ft = (mo/«o) ' is real 
and f & {iufi/gf^ when g y> u>e (the notation is explained in connection 
with (4.24); the g denoting conductivity should not be confused with 
the g denned by (7.20)) the quantity —i xli (/a' s has a positive part. Thus, 
the attenuation of ./„ in the shadow is decreased slightly when the con- 
ductivity g of the cylinder is reduced from infinity to a large finite value. 

7. FIELD AT A GREAT DISTANCE BEHIND THE PARABOLIC CYLINDER 

The field at any point, for the case of horizontal polarization, is 
given b}' expression (5.5) with S->(h) given by (5.4). Since So(h) is the 



* Reference 11, page 424. 
f Reference 5, page 418. 
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only troublesome term in (5.5) most of this section will lie devoted to 
its study. Far behind the cylinder £ and tj are large and positive, and the 
corresponding terms in the integrand of (5.4) are 

T(n + \)U n (zWn(z') = (2it/ v yi 3n e- ir -(l + f)/2 w m , (7.1) 

where the asymptotic expressions (9.16) and (9.17) give 

1 + / - 1 + il ^ r ^ + i(n ± VI" + 2) + <>(»VA (7.2) 
4£- 4?7- 

In writing the "order of" term it is assumed that £ and 77 are both 0(r ) 
with r » n 2 . 

When (7.1) is put in (5.4) we obtain 

SM - V, f ^"^'f,?" «n (7.3) 

J La Bin tt/iTT „Uo) 

where, from the expressions (4.1) for £ and 77, 

M, = [(i/ir) , V ,> sec .(0/2)1/4,, (7.4) 

= Sw/r, = ([tan (5/2)1/7? = cot 0/2 + tt/4) tan (0/2). 

Although it is not proved here, there is good reason to believe that 
(7.3) can be written as 

SM = tf, f ^Zvi\ dn + WW < 7 - 5 > 

Ji t Bin rnW n \Zo) 

when r becomes large and we restrict ourselves to the region | <p \ < ir/2 
in order to get 0(£ 2 ) = 0(?7 2 ) = 0(r). The first term contains r only through 
the factor Mi and is of order r~ m . The "order of" term assumes h to 
be moderately large compared to unity but h 2 « r. When /i < 1 the h 
is to be replaced by unity. 

The general idea leading to (7.5) is that (7.2) may be used over the 
portion of L 2 where the integrand is large and important. On the portion 
where (7.2) differs appreciably from unity the integrand is negligibly 
small. The important portion of L 2 runs from n = % to n — - l A ~ ih 
(approximately). In particular the variation of « 2n V„(^n)/sin m W„(z' ) 
along L 2 may be summarized as follows: from — }^ to + 1 » it decreases 
exponentially as i n , from —Yi to -ih it is equal to — 2i plus an oscillat- 
ing function of order unity, and from — ih to — *« it decreases slowly 
at first and then more rapidly until it goes down like 1 (steepest 
descent behavior). This may be shown with the help of Fig. 12.2, the 
entries for regions I' a and //' in Table 12.3, and the following items [see 
(12.9) and Figs. 10.1 and 10.2 for z = C ll2 no with -3tt/2 g arg (tjg 
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- 2m) < tt/2)]: 

(a) Re \j{t\) — f(to)] is almost zero between — J-£ and — ih. 

(b) Im [/(/i) — f(to)] is almost zero between — ih and —too, 

(c) h/fo runs from zero to unity as n goes from —1 to — 1— ih. 

Items (a) and(b) are consistent with/(* ) - f(h) & (H4)(^ 2 ~ 2m) m 
which holds when n is near — ih and which was mentioned in connection 
with (6.10). 

This concludes our discussion of the reasons for believing that (7.5) 
is true for general values of h. Now we shall check it for the special case 
h = 0. 

When we set h — (i.e. z' = 0) in (7.3), use (9.4) and close L 2 by an 
infinite semicircle, we obtain 



«0)~;-i4r[i- 2 -^+--.J. (7.6) 

This agrees with the first two terms in the asymptotic expansion of 
the Fresnel integral expression (5.20) for S 2 (0). In expanding (5.20) we 
need the first of the two asymptotic expansions (both of which hold 
for T » 1) 

iVa-is^i-^ +...], (, 7 ) 



f T _«i , m i expi-iT^) [ 1 1 

L e *~^/^ + — 2r — L sa* + J* 



and also the first of the relations 



x sin -\- y cos — T% =» — r, 

To = 7,(1 + |8) cos (0/2). 



(7.8) 



(7.9) 



In much of the following work we shall assume £ and t\ to be so great 
that we can neglect the terms denoted by 0(h 3 /r 312 ) in (7.5). We shall use 
the asymptotic sign ~ to acknowledge this omission. 

From (7.4), /3 is equal to unity when <p = 6 — -k/2. When r is very 
large this value of <p marks the shadow boundary. In the shadow /3 > 1 
and in the illuminated region /3 < 1. 

Closing L 2 on the right and on the left converts (7.5) into 

82(h) ~ 2iM 1 £ t3 n V n (zo)/W n (z»). (7.10) 
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&(A) ~ 2iMi j(/3 - I)" 1 - tt 2 [ . - P Z?X /A 1 ) (7.11) 
\ Pt L sm rt dTF„(zo)/dn_]„=„J 

It can be shown that (7.10) converges if < 1 (see (4.18)) and that 
(7.11) converges if > 1 (see (12.13)). The term 1/(0 - 1) in (7.11) 
comes from the poles of esc ini inside the path L 3 shown in Fig. 5.1 
From (7.7) and 

— x sin -4- V cos — T\ = r, 

(7.12) 
2 1 , = 77(1 + 0) cos (0/2) 

it may be shown that when j8 > 1, so that T\ is negative, (5.6) has the 
asymptotic expression 

exp [-ix sin + iy cos 0] + & ~ 2iMi/(l - 0). (7.13) 

When this is added to (7.11) the 1/(1 — /3) terms cancel leaving a series 
for E valid in the shadow where > 1: 

"r ^TntZo) "I , v 

2=f |_sm ttm dW n {zo)/dn_\„ = n, 

This series may also be obtained from the more general series (5.12) 
for E by using (7.1) and neglecting/. 

We now take up the problem of finding the paths of steepest descent 
for the integral in (7.5) when is near unity and h is large. When 0=1 
and h is large, the integrand in (7.5) may be expressed in terms of 
exp \J(ti) — f(k)] by using Table 12.3. In Section 6 it has been pointed out 
that the path of steepest descent for exp \f(h) - f(U)] is the path ACD 
of Fig. 6.2, with C being the high point. This suggests that the path 
ACD should be used to deal with the terms in (7.5) leading to exp [/(h) 
- /(ft,)]. These terms are U n (z' )/W n (z' ) (introduced by the use of (4.8)) 
for the portion of L 2 between B and C, and V n (z' )/W n (z' ) for the portion 
between C and D. As a further argument supporting the use of the path 
AC we note that when n is on AC, i.e., on the edge of region I'b, Table 
12.3 gives 

U n (z' )/W n (z' ) ~ -i(l - i^Xh/to) 1 ' 2 exp \f(h) - f(t,)) . (7.15) 

Hence the variation of i n esc im in the integrand of (7.5) is just cancelled 
by that of (1 - i -4 ") in (7.15). Consequently i 2n U„(z' ) /[sin *nW n (z' Q )] 
varies as exp [f(ti) — /(M] along AC (the variation of /i// is relatively 
small). 

These considerations lead us to write (5.4) as 



462 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1954 

S- 2 (h) = Sn + S22 + &3 
S 21 = - [ F dn, 

J B 

S T2 = -J FU„U) dn/Wniz!,), (7.16) 

S 2 , = [ FV n (z' ) dn/W n (z' ), 
Jc 

F = e*sec (6/2)(iw/2)"T(?i + l)U n (z)W n (z')/2i sin rn. 
When instead of (5.4) the expression (7.5) for 82(h) is used we obtain 

S21 ~ —Mi I i 2 "0" dn/ sin mi, 

J B 

Sn Mi f i 2n p n U n (z' ) dn/[sm wnW n (z' )], (7.17) 

Sn ~ Mi / i 2n p n V„(z' ) dn/[sin irnW n (zo)l 
Jc 

In $22 it is permissible to swing AC from its original position BC 
because the zeros of U n (z' ) cancel those of sin xn. When |3 = 1, AC and 
CD are the paths of steepest descent for £22 and S23 in (7.17) because 

Im lf(h) ~ /(Ml - on ACD. 

The asymptotic expression (7.17) for Sn may be evaluated by tem- 
porarily assuming to be a complex number with | | < 1 and | arg | 
< 7i72. The integral along BC is the integral along BCE minus the in- 
tegral along CE (see Fig. 6.2). Deforming BCE into L x of Fig. 5.1 shows 
that its contribution to <S 2 i is — 2iMi/(l — 0). An infinite series for the 
integral along CE may be obtained by expanding i "/sin irn in powers 
of exp ( — i2irri) and integrating termwise from n = n = —1 — ih to 
n = 00 — ih, i.e., from C to E. In this way we obtain 



&i ~ -2iMi 



1 ■A exp (no log — 27r//i) 



(7.18) 



_1 - 1=0 log - 2in7 

Despite the appearance of the right hand side, it is analytic around 
= 1 and analytic continuation may be used to show that (7.18) 
holds for < < w. 

When h is large only the first term in the series is important and 
we have 

Sn ~ 2iAfi[(0 - I) -1 - /r 1_, 7log 0] 

(7.19) 
= 2iM 1 (fi - l)" 1 + iM 2 g x 
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where we have introduced two quantities which will be used later: 



0" h - 



exp [—ir + igh ll3 \ 



*.-»*V-- W 2€ sin(,/2) • (7.20) 

g = -/i 1/3 log/3. 

When |8 = 1 

8 Z1 ~ 2iMi(ih'+ }4). (7.21) 

When h is large most of the contributions to the integrals (7.17) for 
£22 and S 2 3 come from around n = n = —1—ih, and we may use the 
approximations 

U n (z )/WM ~ C^AHaT^/AHa), (7.22) 

F*&)/TP»«) ~ FAiiat^/Aiia), (7.23) 

a = (ihy m (n + 1 + ft), n - Wo - a(ih) m , 

which come from (13.21). Setting these in the integrals (7.17) and using 
the fact that t 2 "/sin irn is nearly 2i around n leads to 

Sn M 2 f exp( -i m ag)Ai(ai~ m ) da/Ai(a), (7.24) 

Jooexp (»2jt/3) 

-Mcxp <— 12»/3) 

£ 23 - -t 2/8 M 2 / exp (-i V3 ag)Ai(ai 413 ) dot/Alia), (7.25) 

Jo 

& + S n ~ iMMg), (7-26) 

where 

*( ff ) = i [ exp (-i llz ag)Ai(.m^ 13 ) da/Ai(a) 

Jooi4/3 

+ i 5/3 f exp (-i m ag)Ai(ai 4,i ) da/Ai(a). 

Jo 

The expression (2.11) for V(g) is obtained from (7.27) by changing the 
variables of integration and using the transformations (13.17) for Ai(a). 
Thus, when h is large and /3 near unity, (7.19) and (7.2G) give 

St(h) ~ 2il\h((S - l)" 1 + iMJg~ l + ¥(?)]. (7.28) 

In the shadow, where (3 > 1 and g is negative, (7.13) and (7.28) give 

E = exp[-u: sin + iy cos 0] + S, + S,(h) (7.29) 

^iaf.to -1 4-^)]. 

This and the series (7.14) for E suggest that ¥(g) 4- \/g may be ex- 
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pressed as a series in which the parabolic cylinder functions in (7.14) 
are replaced by Airy integrals. One way of obtaining this series from 
(7.14) is to use the Airy integral approximations (13.21). The zeros 
rii, n 2 , • ■ • of W n (z' ) go into the zeros a h a 2 , • ■ • of Ai(a) by virtue of the 
relation n — n = a(ih) and we have 

E~ l ^M i ± eX v}~ a Z :S] , (7.30) 



[Ai'{a s )Y 



.-l/3i 



■9(g) + 1/g = -i 2^ rjTrTrTii » t'- 31 ) 

s =i [At {a s )\ z 

where g < 0. Here, as in (6.17), a x = -2.338 • • and Ai'{a x ) = .701 • • • 
In obtaining these relations we have used t 2n /sin irn tt 2i, and 

Ai(a/") = -i hlz Bi{a,)/2 = i 5 ' 3 /2tAi' (a.) , (7.32) 

where the first equation follows from (13.17) and the second from the 
Wronskian 

Ai(a)Bi'(a) - Ai'{a)Bi{a) = 1/tt. (7.33) 

The equal sign in (7.31) holds even though the steps leading to it 
indicate that ~ should be used. This may be seen from an alternative 
derivation of (7.31) in which Ai(oa ) in the first integral of (7.27) is 
replaced by the right hand side of* 

Ai(oT 4 ' 8 ) = -i A,3 Ai(a) - i-^Aiia*?' 3 ). (7.34) 

In the first portion the Ai(a)'s cancel and the resulting integral con- 
tributes — 1/g to (7.27) (g must be negative for convergence). The 
second portion combines with the second integral in (7.27) to give a 
contour integral which leads to the series in (7.31) when the path of 
integration in the a-plane is closed on the left. The closure may be 
justified by the asymptotic expressions (13.19) and (13.20) for Ai(a) 
(again g must be negative). 

Since the integrals in (7.27), and their equivalents in (2.11), converge 
uniformly for all finite values of g, V(g) is an integral function of g. 
When g is negative ~^(g) may be computed from the series (7.31). When 
g is positive I have not been able to find a practicable method of ob- 
taining ^((7) other than the numerical integration of (2.11). The results 
are shown in Table 2.1. Since ~&(g) is an integral function its Taylor's 
series about, say, g = — .5 converges for all values of g. The coefficients 
in this series may be computed from (7.31). However, I was unable to 

* Reference 11, page 424. 
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obtain useful results by this method because the computation of the 
coefficients becomes more and more difficult. 

When g is large and positive it may be shown that 

*(</) ~ - f l + (irg) m exp (-tf/12). ( 7 - 35 ) 

The procedure used to establish (7.35) is much the same as that used 
to establish the more general result 

£22 + £ 23 ~ -iMog~ l 



7i(l - 0)' 



1/2 



2 exp[-w- + 2ih(l - j8)/(l +0)] 



_r(l + j8)_ sin \{<p + B + tt/2) ' (7.36) 

1 - g = sin |(g - 9 + tt/2) 
1 + sin |(v + + t/2) ' 

which holds when /i 1/3 (l — /3) is large and positive. 

When <p is near —tt/2 + 0, (7.36) gives the same result as (7.26) 
plus (7.35). For <p near -tt/2 + d, 

g & h m (l - 0} « [2/i 1/3 /sin 0] sin i L - -f ^j , (7.37) 

which shows that g is proportional to the cube root of the radius of 
curvature 2/t/sin 3 9 at the point where the incident ray is tangent to the 
cylinder. 

When j8 < 1, (7.36) may be obtained from 

£•■•> + &. Mt / -4-^- dn - Mi / ^7-TK dn. (7.38 

J c sin 7m J.4CB sin 7mlr„Uo) 

The second integral in (7.38) represents, asymptotically, the wave 
reflected by the cylinder. This interpretation is suggested by the fact 
that, when the expression (7.1) for T(n + l)U n (z)W n (z') is substituted 
in expression (5.1) for E, the resulting integral may be written as the 
second integral in (7.38). 

The first term in (7.36) is obtained when i 2n /sin irn in the first integral 
of (7.38) is approximated by 2i and the result integrated. When the in- 
tegrand of the second integral is examined with the help of Table 12.3, 
it is found to have a saddle point* at m = mi on the imaginary axis 
between m = and m = —ih. Near w?i the integrand is approximately 

2/T 1 (/ 1 //o) 1/ ' 2 exp [F(m)] (7.39) 

* It is interesting to note that a saddle point also appears in the study of re- 
flection from a sphere. See page 86 of reference. 7 
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where F(m) = f(h) - f(t ) + to log p and F'(to) = log ikP/k). Here 
/ and h are functions of to defined by (12.9). At TOi we have fo/3 = *i 
and this leads to 

Tin = -4^/(1 + /3) 2 , F(na) = 2ih(l - p)/( 1 + fi), 

F"( mi ) = -(1 + 0)M(1 - /3). (7.40) 

When we attempt to deform the path of integration ACE of the 
second integral in (7.38) into a path of steepest descent, we encounter 
no trouble near TOi in regions I' a and I'b. The path passes through Wi 
with arg (dm) = — t/4. Soon after passing through to-i the path of 
steepest descent strikes the boundary between I' a and //' at a point 
we shall call G. At G the imaginary part of to is 2h(l — &)/(l + fi) log p. 
The asymptotic approximation to the integrand changes its form at 
this point. The choice of the path from G out to °o is not important 
since it contributes little to the value of the integral. However, if we 
insist on following paths of steepest descent, it turns out that we must 
split the path of integration at G. 

When h m (l — P) » 1, it may be shown that the value of second 
integral in (7.38) is nearly equal to 

-2M 1 [-2Tr/pF"(m 1 )} 1 ' 2 exp [F(toO] 

and this gives the second term in (7.36). 

So far, in this section, we have been dealing with the case of horizontal 
polarization. Since the work for the case of vertical polarization (in 
which H plays the role of- the wave function) follows much the same lines, 
we shall merely list the formulas corresponding to those already ob- 
tained for horizontal polarization. M h p and M 2 , g are still given by 
(7.4) and (7.20), respectively. 

sm = m, f \ - 7£rr\ dn + °^ 3 A 3/2 ), ( 7 - 41 ) 

Jl-, sm irnW n {zo) 
S 3 (0) = -&(0), (7.42) 

SM ~ 2iM 1 Z P" 'V n (z'o)/'W n (zo), P < 1 (7.43) 



S 3 (h) ~ 2iMi <(P - 1) '-iL a<TT7 ; /x /a 

' [ f=\ sin wn d W n (zo)/dn_ 



P > 1 (shadow region), 
n» — sth zero of 'W n (za), 



(7.44) 
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H ~ -2iMnr E 



■ r"/3" 'VJjzj) 
sin 7m 3 



^#^1 , 0>1 (7.45) 

Trn(zo)/dn_| n = n ; 



SxQl) = $31 + $32 + $33 , 

$31 = $21 defined by (7.16), 

S, 2 = $,, with 'Unizid/'Wniz'o) in place of U n (z' )/W n (z' ), 

$ 33 = $23 with 'V n (zo)/'W n (zo) in place of V n {z'o)/W n {z',), 



(7.46) 



* 2/3 M 2 /"' 



$32 ~ C % / exp (-t^Ui'CoT 4 ") da/Ai'(a), (7.47) 

>oo exp («2ir/3) 



.oo exp (-i"2ir/3) 



$33 ~ M 2 / exp {-P 9 ag)A?{ar) da/ At' (a), 

Jo 



$32 4- $33 ~ il\m v (g), 



(7.48) 
(7.49) 

(7.50) 

(7.51) 

g-^t r^yZ ^<o (7.52) 



tf.fo) = r 1/3 f exp (-i 1/ 'a fl r)Ai # (ar 4/ ') da/ At' (a) 

Jooexp (i'2)r/3) 

»oo exp (— 12*73) 

-i / exp (-i 1/3 a^)^t 7 (a? 4/3 ) da/Ai'(a) t 

Jo 

H ~ Mf.fo"" 1 + * v (g)], g < 



fe (-a^AtCa!)] 



' .-l/3\ 



*.M + g- 1 = -^ 1/3 Z 



/ ' i/a\ 

exp(-a.flri ) 



< (7.53) 



h(-a a )[Ai(a s )] 2 ' 
a. = sth zero of Ai'(a), a[ = -1.019, Ai(a[) - 0.5357, 
At'CaiO = -i in Bi'{a,)/2 = -i""/©^^)], (7-54) 

Ai"(a) = a/h'(a). 

When g is large and positive, 

%(g) ~ -g~ l - d*g) ui ™p (-tf/M), (7.55) 

$32 + $33 ~ -iMig~ x 

"h{\ - /3) 



+ 



1/L> exp [-ir + 2*(1 - |8)/(1 + 0)] (7.56) 



_r(l + fi)_ sin |(<P + + r/2) 



The change in sign of the second term on the right in going from (7.36) 
to (7.56) comes from (12.2) and the analogous expression for 'W n (z' ) 
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(only h contributes to U n (z' ) and only to to W n (z' ) at the saddle point 
mi of the second integral in (7.38)). 

So far in this section the parabolic cylinder has been assumed to possess 
infinite conductivity. When the cylinder has a finite (but very large) 
conductivity, it may be shown that the field far out in the shadow is 
approximately 

r, ,r f ***" V^ + ^'Vnti) , 

E ~ Mi \ T „ , u — -i nxr ( /v dn. (7.57) 

J Li sin irn W„(z ) + a 'WJ&l) 

Equation (7.57) is suggested by (7.14) and (4.25). The analogue of 
(7.57) for vertical polarization may be obtained by replacing E, a in 
(7.57) by H, t so that 'V„(z' ) + TV n (z' ) appears in place of V n (z' ) 
+ <r~ 'V(z f ), and so on. 

When the parabolic cylinder functions are replaced by Airy integrals 
according to (13.21) and (13.24), equation (7.57) may be written as 

m„ f [exp {-agi m ))Ai[{a 4- k)i 4 '*] 



E ~ ("'Mi 



f [exp (-agin\Ai[( a + kK'l A (ej Ksn 

U Ai(a + k) ^ (? - 58) 



where g and M 2 are given by (7.20), a by (7.23) and 

h = -(ih)~ m $/U. (7.59) 

| k | is small compared to unity. L 4 is a path of integration in the a plane 
which encloses the zeros of Ai(a 4- fc) in a clockwise direction. Changing 
the variable of integration in (7.58) to u = a + k enables us to con- 
clude that 

I" E for finite 11" / _t£\\ \E for infinite] , ? . 

[_conductivityJ ™ |_ \ fo / _ _conductivity J ' 

Since we have assumed 5 = t/2, the relation (7.60) holds in the region 
where the angle \f/ defined by Fig. 2.3 is negative. 

The analogue of (7.58) for vertical polarization is obtained by re- 
placing E by H, omitting the i m , and replacing the ratio of the Airy 
integrals by 

Ai'[(a + t/a)?'*\/Ai'(cL + 1/ a) (7.61) 



where 



vl/3 



1= -(ihy'Wta. (7.62) 

Even though h is large, f/fo is assumed to be so small that t is small 
compared to unity. The path of integration L 4 must now enclose the 
zeros of Ai'(a + I fa) which are close to those of Ai'(a) at a = a' s ,s = 1, 
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2, • • • .It must not pass close to a = since the work leading to (7.61) 
assumes (fa to be a small number. Changing the variable of integration 
to v = a + C/a, approximating (/a, {/a by f/v, f/v 2 , and evaluating 
the integral by considering the residues of the poles at v = a, gives 

h M r**, ±(i- -Cf ' exp 'r'"y;-,^f " :)l • (7.63) 

•-i \ a,"/ { — a s )[Ai{a s )\~ 

This shows, to a first approximation, how the expression (7.52) is 
modified when the cylinder is a very good, but not perfect, conductor. 
Of course g must be negative in (7.63). Since ( in (7.62) varies as li 3 
while k in (7.59) varies as h~ in it appears that the field for vertical 
polarization is much more sensitive to changes in the conductivity 
than it is for horizontal polarization. 

It may be verified that the change in the exponential terms in the 
series (7.30) and (7.52) produced by finite conductivity, namely 

a s changes to a 8 — k 

' i. t ' « ' (7 ' 64) 

a s changes to a, — €/a„ 

agrees, to a first approximation, with the change produced in the cor- 
responding series (given, for example, by the series (27) and (28) on 
page 45 of Reference 7) for the propagation of radio waves over the 
earth's surface. 

8. FIELD AT A GREAT DISTANCE BEHIND THE PARABOLIC CYLINDER WHEN 
6 = 7r/2 AND h IS LARGE 

In the work of Section 7 the angle of incidence 8 may he anywhere 
between and ir. Here we take 6 = tt/2, which corresponds to the case 
shown in Fig. 2.3 and described in Section 2. Some simplification is ob- 
tained thereby. For example, the incident wave is now simply exp ( — ix). 
We shall write the expressions for the horizontal and vertical polariza- 
tion cases as 

E = (e"*" + &)r + Sa + (£22 + &■), (8-1) 

H = (e- ix + Si), + S n + (£32 + £„), (8.2) 

respectively. Here S 2h • • • are defined by (7.16) and (7.46) in which 

6 = tt/2, iv = 1, 

(8 3) 
= Z/v = cot {tp/2 + tt/4) = 1 - <p + <p/2 - v/S + • • • 

Throughout this section will be defined by (8.3), i.e., by (7.4) with 
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6 = jt/2. Also, from (5.3) and (5.6) 

(e~ iz + ft), = (t/rW* r e- u * dt, 7\ = 2~ l % - |) 

■*-- (8.4) 

= (2r) 1/2 sin Gp/2). 

The subscript r is used to denote correspondence to a half-plane with 
its edge at r = 0. 

When h is large, physical reasons lead us to expect a similarity be- 
tween our field and the one behind a half-plane with its edge at the 
crest of the cylinder where p = (see Fig. 2.3). The main part of this 
field is the analogue of (8.4): 

(«-*■ + ft), = W*) W ** \2 *** di, T 3 = (2p) 1/2 sin (*/2), (8.5) 

jhere the subscript p indicates that [exp (—ix) + *Si] p corresponds to 
Jiffraction behind the half-plane just mentioned. 

In order to make use of the similarity between the field behind the 
cylinder and the half-plane with its edge at the crest of the cylinder, 
we change the polar coordinates from (r, <p) to (p, \p). From 

pe 4 * = re'* - ih (8.6) 

it may be shown that, when h 2 /r is small, 

(e~ ix + &)r - («~* + ft), - (t'AO'V* f V" 

[e'"" f -1] r 0(/i 3 /r 3/2 ) 



(8.7) 
= 2iMi i pi h sin v _n _,_ nft.*/..»i*\ 

where Mi is obtained by putting 8 = 7r/2 in (7.4). 

When we combine (8.7) and the expression (7.19) for 5 2 i the 
2iMi/(fi — 1) terms cancel leaving 

(e~ ix + Si) r + Sn = («"* 4- ft), + |^- e rt 8in ' + i^ 2 

0-1 (8.8) 

+ 0(h'/r 312 ) + 0[Mi exp (-2ich)]. 

The sum of the terms involving M\ and Mi may be expressed in a 
form which contains the expression c(r) defined by (2.5) and the quan- 
tity b defined by 
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b = -log = log tan (J + |) = 9 + *'/6 + • ■ • 

,-i/3„ (8.9) 

tanh b = sin <p. 

Replacing c{r) exp (i/i sin ^j) by c(p) plus a correction term then con- 
verts (8.8) into 

{e- U + Si) r + Sn = («"*• + S X ) P 

, c( P )[l + exp (26,)] I/2 |~_L exp gfeb - gt tanh 6) ~| , m \ 

+ yir- - Lr^7 + 1 J p v- 10 ) 

+ 0(ft 8 /r 3/2 ) + Olr-V^, 

where the subscript p on the square brackets indicates that b is to be 
replaced by b p denned by 

b p = log tan (*/2 + tt/4) = <A + * 8 /6 4- • • • (8.11) 

The quantity within the square brackets in (8.10) is continuous at 
6 = where it behaves like (neglecting 0(6) terms but retaining 0(hb 2 )) 

Y 2 4- ihtf/Z = l A + ih ll3 g 2 /3. (8-12) 

Expression (8.10) is to be used with (£22 + £23) and (£32 + £33) ob- 
tained from (7.17) and (7.46) (with = ir/2 and w = 1), respectively. 
When \f/ is small, expression (8.10) becomes 



(e~ ix + S 1 ) r + Sn = (e~ ix + &), 
+ c(p) 



1 _ 1 h m exp (V/3) 1 ( 8 - 13 ) 

2 ^ + <? J P + 



The subscript p on the square bracket indicates that g is to be replaced 
by g p defined by 

g , = h ll3 b p = h m (+ + * 3 /6 +•••)• (8.13) 

When, in accordance with (8.1) and (8.2), we add to (8.13) the ap- 
proximations (7.26) and (7.49), namely 

S n + & 3 ~ *M&(g), 

(8.14) 
£32 + S 3Z ~ iM&M, 
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we obtain 



E = (e~ ix + Si), + c(p) fi - i + <f- + *(g)\h m exp (i//3)] 



+ 



H = («-*■ + ft), + c(p) 



+ <- - +*,fo)K exp(V/3) 



(8.15) 



(8.16) 



+ 



The terms neglected in (8.15) and (8.16) are the "order of" terms in 
(8.10), plus those neglected by virtue of $ being small, plus the errors in 
(8.14) The errors in (8.14) are of two kinds namely those of 0(/i 3 /r 3 2 ) 
and those due to approximating the parabolic cylinder functions by 
Airy integrals. 

It is interesting to observe the forms assumed by (8.15) and (8.16) 
when h = even though they are not supposed to hold for small values 
of h. In this case p, \f/ go into r, <p and the right hand sides of (8.15) and 
(8.16) become the same, namely 

(e~ ix + Si) r + c(r)/2. (8.17) 

The half plane results given in Section 2 become, for small values of <f, 

|| = (a - * + SO, ± c(r)/2 (8.18) 

where the upper sign corresponds to E and the lower one to H. Com- 
parison of (8.17) and (8.18) shows that (8.15) for E reduces to the proper 
value but (8.16) for H fails to do so because the signs of c(r)/2 do not 
agree. 

The discrepancy is apparently related to the approximations we have 
made in obtaining the expression (7.19) for £ 2 i from (7.18) and to the 
errors introduced by approximating the parabolic cylinder functions 
by the Airy integrals. As we let h — > in the more complete expression 
(7.18) for $21, the value obtained for £ 2 i — > °° . This is explained by the 
fact that the upper limit of integration — 1 — ih (at point C) approaches 
the pole of the integrand of (7. 17) at n = — 1 . This large value of Sii tends 
to be cancelled by the large value of £23 (for horizontal polarization). 
On the other hand our approximation (7.19) yields via (7.21) the value 
iM\ for £21 and S 3i when h = and = 1. The factor h m in il/ 2 makes 
our approximations for £22, £23, S32, £33 in terms of Airy integrals vanish 
when h = 0. 
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Incidentally, if, instead of taking the point C of Fig. 6.2 to be at 
— 1 — ih, we take it to be at — 3^ — ih (a choice which receives some 
support from the Airy integral representation obtained from the view- 
point of the differential equations discussed in the first part of Section 
13), the approximate integrals of (7.17) and (7.46) may be integrated 
directly when h = and /3 = 1. It is found that 

8n~~— log 2, S* ~ ** log 2 + **1 , 8» ~ +iMx/2, 

IT IT 2 

S 81 ~ _ — log 2, £ 32 ~ — log 2 - — - , S 33 iMi/2, 

T 7T 2 

and these add to give the values &(0) ^ iftfi, #3(0) ~ —iM\ required 
by the half- plane case. 

It is seen that a rather thorough investigation of the errors introduced 
by our approximations would be required to resolve the discrepancy 
between (8.17) and (8.18). Since we do not intend to go into this subject, 
and since the errors we have made may be as large as the % which ap- 
pears within the square brackets of (8.15) and (8.16), we shall "split 
the difference" between the two polarizations and omit the Yi. alto- 
gether. This is done in Section 2 where t ~ g p . 

9. THE FUNCTIONS U n (z), V„(z), W n (z) 

The functions U„(z), etc., are defined for all values of z and n by the 
integrals 

u n (z) = J-. f r"-V-' 2+ - di, 

2m Ju 

Vn(z) - J-. ( r-V" 4 * 1 dl, (9.1) 

Zirl Jy 

W n {z) =-L. f r^e- Ml dl. 

2m Jw 

where the paths of integration U, V, W in the complex {-plane are shown 
in Fig. 9.1. The cut in the {-plane runs from — =0 to and has been 
introduced in order to make the function t~"~ l one-valued. In some of 
the later work the paths of integration will cross this cut. Of course, 
this requires close attention to arg t. 

The initial and final points of the various paths (denoted in Fig. 9.1 
by the subscripts i and /) are located at infinity. Arg t = -ir at (/,■ 
and W f and + t at U f and F,-. 

We shall give a summary of the properties of the functions (9.1) 
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which will be needed in our work. These functions are related to the 
parabolic cylinder function D n (z) ls ' 19 through the equations 

U n (z) = 2" /2 e'" 2 D n (2 m z)/T(n 4- 1), 

VM = -t- r 12 e* V2 D_ n _x {-iz2 m )/{2^) m , (9.2) 

W n (z) = -i n 2 nl2 e' 2 ' 2 D. n ^{iz2 m )/{2^ l \ 

We use the functions U n (z), etc., here instead of D n (z) because they 
seem to be more convenient for the particular problem we have to 
deal with. 

From the definitions (9.1) it follows that U n (z), V„(z), W n (z) are one- 
valued analytic functions of z and n. By expanding exp (2zt) in (9.1) and 
integrating termwise it may be shown that 

U n (z) = 2 A cos (irn/2) + 4=zB sin (tu/2), 

V n (z) = -Ar n -2zi~ n+1 B, (9.3) 

W n (z) = -Ai n -2zi n ~ x B, 

When z = and U' n (z) = d U n (z)/dz, etc., 



r(l + n/2) ' 2r(l+n/2)' 

Wn(0) = 



— I 



2r(l + n/2) ' 
v , m 2 sin (W2) r m -T" , 0.4) 

W' n {0) = 



W 



18 See E. T. Whittaker and G. X. Watson. Modern Analysis, Fourth Edition 
(1927) Cambridge Univ. Press pp. 347-354. 

19 W. Magnus and F. Oberhettinger, Formeln und Satze ftir Speziellen Funk- 
tionen, 2nd Ed., Springer, 1948 Chap. 6 Section 3, and p. 227. A comprehensive 
account of D n (z) is given in the forthcoming work, Higher Transcendental Func- 
tions, compiled by the staff of the Bateman Manuscript Project. 
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Let T n (z) denote any one of U n (z), V n (z), W n (z), let primes denote 
differentiation with respect to z, and let asterisks denote complex con- 
jugates. Then we have the following relations 



T"(.z) - 2zT' n {z) + 2nT n {z) = 0, 

£ tr'T'M] + 2ne-''T n (z) = 0, 

az 

f, [e-'" 2 T n (z)] + (2n + 1 - f)*-* n TM - 0, 

dz- 

T' n (z) = 2T n - X (z), 

7iT n (z) = 2zT n . 1 (z) - 2T n ^{z), 

U' n (z) V n (z) - U n (z) V'M = «'2"e-' 2 A 1/2 T(n + 1), 

71(f) W n (z) - V n (z) W' n {z) = i2 n e z *fr m T(n + 1), 

W' n {z) U n (z) - W n (z) U' n (z) = i2V 2 A 1/2 T(n + 1), 

U n (z) + V n {z) + W n (z) = 0, 
[V n (z)V = WAz*), [W n (z)\* = F..(«*), 
[U n (z)]* = UA**), 
V n (-z) = T 2 " W n (z), W n {-z) = i 2n V n (z), 
U n (-z) = -i n V n (z) -%** W n (z), 



V-n-tiiz) = -i" +l KU n (z), 



(9.5) 
(9.6) 

(9.7) 
(9.8) 

(9.9) 

(9.10) 

(9.11) 

(9.12) 




Fig. 9.1 — Paths of integration used in the integrals of equations (9.1) which 
define the functions U»(z), V„(z) and W n (z). The subscripts i and / stand for the 
"initial" and "final" points of the paths. 
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U-^iiz) = f*V - i- 2n )KW n (z) = -(2t)-V /2 e -'V n ( Z )/r(-n), 

U- n -,(-iz) = i n ~V n - C 2n )KV n {z), 

K = e- zl Y(n + l)/x 1/2 2 n+1 

Equation (9.5) may be obtained from (9.1) by forming T*(z) — 2zT'(z) 
and integrating by parts. Equation (9.10) follows from (9.1) upon 
joining the paths U, V, W to obtain a closed path of integration. From 
(9.11) it follows that when we have an expression for V n (z) which holds 
for all values of z and n, replacing i by —i (or f ) gives the correspond- 
ing expression for W n (z). Equations (9.13) may be obtained by using 
the fact that U- n -\{iz) exp {z") etc. are solutions of the differential 
equations (9.5) and (9.7). 

The relations (9.11) and (9.13) enable us to compute the values of 
U n (z), V n (z), W n (z) for z = i in p and z = i~ 12 p and all n when the 
values of any two are given for z = i m p and Re(n) ^ — %. 

It may be verified that as z becomes large and n remains fixed the 
differential equation (9.5) has the asymptotic solutions 



. , 2V „ / n 1 - n , . a 



-n-1 z* 

2 e 



S2 ^ w ' ^ = 2\Ar 2i? ° ' ° ' ° ' ;1 ^ 2 ' 



(9.14) 



where 



| arg z | < t and s 2 (— n— l,iz) = i n Ksi(n,z), 

K being given by (9.13). In terms of these functions we have 

V n {z) rsj —is 2 (n, z), < arg z < ir 

~ — si(n, z) —isi(n, z), —tt/2 < arg z < (9.15) 

~ — Si(n, z) — i~ 4n+1 s 2 (n, z), — -k < arg z < —tt/2 

PT„(z) ~ is 2 (n, z), —k < arg z < (9.16) 

tfn(z) ~ si(n, z), -tt/2 < arg z < tt/2. (9.17) 

The first expression for V n (z) in (9.15) follows when we note that the 
leading term may be obtained from (9.1) by choosing the path of in- 
tegration V to be t = z + t where t runs from -co to + oo t and | z \ 
is supposed to be large. (9.17) follows from the first of (9.15) and the 
relation (9.13) between V- n -i(iz) and U n (z). Asymptotic expressions 
for W„(z) may be obtained by taking the conjugate complex of those 
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for V„{z). The second and third expressions for V n {z) follow from the 
other asymptotic expressions and (9.11), (9.12). 

When R(ri) < 0, the theory of gamma functions and (9.1) lead to 

T(n + l)r(-n) U n (z) = -* esc irriU n (z) 

r 80 (9.18) 

= / r" n_1 exp (-• t 2 -2tz) dr. 
Jo 

By expressing v^ exp [-(z — tf] as the integral of 

exp[-r 2 + 2i(z - t)r] 

taken from t = — oo to + °° and substituting in (9.1) it may be shown 
that, when R (n) > —1, 

U n (z) =Fi n r e- t2 - 2Ut t n dt, 

J— QO 

n ( 2 ) = -Ft— C e- r2+2i " T n dr, (9.19) 

Jo 

(z) = -Fi" f e' r2 - 2i " r n dr, 
Jo 



V 

w 



F = 2V/r(n + l)7r 1/2 . 



When n is not an integer the path of integration in the integral (9.19) 
for U n (z) is indented downward at the origin. Equations (9.19) mav 
also be obtained from (9.1) by using (9.13) and (9.18). 
When n is an integer 

Un (-z) = (-)" U n (z), V n (-z) = (-)" W n (z), (9.20) 

and when n is a positive integer 

u n (z) =sMz) = (-y c -f n e-* 2 , 

ft! dz n 

U- n (z) = 0, (9.21) 

VUz) = -WU*) = -is,(-n,z) = --±- n ^,\ 

From Maclaurin's expansion and (9.21), 

£ t n U n {z) = exp [-t 2 + 2zt\. (9.22) 

n=0 
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10. FORMULAS FOR THE SADDLE-POINT METHOD 

Much of our work involves the behavior of the parabolic cylinder 
functions as functions of n when n is a large complex number. Although 
this subject has been studied by several writers, 20, 1 " 22 their results 
are not in the form we require. As the work of Sections 6 and 7 shows, 
the paths of steepest descent for the integrals in our electromagnetic 
problem are intimately connected with the function f(t ) — f(ti). In 
turn, this function is closely related to the saddle point method of evalu- 
ating U n (z), etc., for large values of n. For the sake of completeness, 
we shall outline this method. We "shall pay special attention to the rela- 
tive importance of the two saddle points as n moves about in its complex 
plane. 

When we write the integrand of the integrals (9.1) as exp [f(t)] we 
obtain expressions of the form 



UM = i f exp [/(/)] dt, 



(10.1) 
f(t) = -I 2 + 2zt - m log i, to = n + 1. 

The saddle points of the integrand are at the points U and t\ in the com- 
plex t -plane where f'(t) is zero: 

2fo - 2zto + to = 0, to - z( = -to/2, 



. z + (z 2 - 2m) m ... 

'o = s ' to + h = Z, 

t z-(z 2 - 2m) m 0i . 

h = 1 — ^ — , Mod = to. 



(10.2) 



Let the path of integration U of (10.1), for example, be deformed so 
as to pass through a saddle point, say t , along a path of steepest descent. 
Let 

fit) = /(/„) - J b,(t - Uf/kl. (10.3) 

2 

Then, if 6 2 is not too small, the contribution of the region around to 

20 Nathan Schwid, The Asymptotic Forms of the Hermite and Weber Functions, 
Amer. Math. Soc. Trans. 37, pp. 339-362, 1935. References to earlier work will be 
found in this paper. Schwid's work is based on R. Langer's study of the asymptotic 
solutions of second order differential equations. 

21 O. E. H. Rydbeck, The Propagation of Radio Waves, Trans, of Chalmers Univ. 
of Tech. 34, 1944. 

22 G. N. Watson, Harmonic Functions Associated with Parabolic Cylinder 
Functions, Proc. London Math. Soc. (2) 17, pp. 116-148, 1918. 
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to the value of the integral is exp [/(* )] times 

dt 
~ (2t6 2 )" 1/2 [1 + j -b£, + lObfo} (10.4) 

-6«£ 3 + [35b| + 566 3 6 6 ]-B 4 - 21006l6 4 B B + 55(280)foiB 6 } 



^/«qi[-Z6t(<-«D)V*> 



where B k = (2b 2 )~ k /k!. The sign of (2ir6 2 )~ 1/2 is chosen so that the argu- 
ment of the right hand side of (10.4) is equal to arg (dt) at t = to on 
the path of steepest descent. The derivatives of f(t) at tg give 



b 2 = 2(' - /i)// , b 3 = 4/i//o, b, = -12h/tl 



-b&2 + mlB 3 = 



h(ti + 9'p) 
24t (t Q - tt)* 



(10.5) 



The values of these quantities at the saddle point h may be obtained by 
interchanging * and k. If more terms of (10.4) are desired they may be 
obtained from the formal result 



o~ / exp 

ZTT J-00 



- !>,//*;! 



dt 



(2ira 2 ) 



-1/2 



_ (10.6) 



1 + £ F»(0, 0, -a ;i , -04, • • • , -a**)/* !(&*)' 

A- =2 



where F n (a tj a 2 , . . . , a„) is the Bell exponential polynomial. 23 It is neces- 
sary to rearrange the terms given by (10.6) in order to get them in 
groups having the same order of magnitude. A more careful treatment 
of the terms in the asymptotic expansions for D„(z) has been given by 
Watson. 22 His method is similar to that used by Debye for Bessel func- 
tions. 

In our work we shall deal with two different complex planes, and the 
reader is cautioned against confusing them. One is the complex i-plane, 
shown in Fig. 10.1, which contains the paths of integration for integrals 
such as (10.1). The other is the complex m-plane, shown in Fig. 10.2, 
which is introduced because we are often more interested in U n (z), etc., 
as functions of m = n -4- 1 than as functions of z. In the earlier sections 



23 E. T. Boll, Exponential Polynomials, Ann. of Math. 35, pp. 258-279, 1934. 
The polynomials are tabulated up ton = 8 by John Riordan, Inversion Formulas 
in Normal Variable Mapping, Annals of Math. Stat. 20, pp. 417-425, 1949. 
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we have spoken of the complex w-plane, but this is essentially the 
m-plane shifted by unity. 



Since we are going to deal with a fixed value of z (i £ or i~ j/) but 
with a variable value of m, we make to and h one-valued functions of m 
by cutting the m-plane as shown in Fig. 10.2. 

It may be shown that to and t\ lie in the opposite half-planes in- 
dicated in Fig. 10.1. This restricts arg U to lie between arg z — 7r/2 and 
arg z + 7r/2. Arg h is restricted to lie between arg z — t and arg z + t 
by the cut shown in Fig. 10.1. It may also be shown that 



k\ ^ Ul 



arg £ — arg 1 1 1 ^ ir. 



(10.7) 



t REGION 




BOUNDARY BETWEEN 

t HALF-PLANE AND 

t, HALF-PLANE 

Fig. 10.1 — Diagram showing the half -plane regions to which the saddle points 
to and <i are confined in the /-plane. 

One might wonder why cuts in the m-plane are required since it has 
already been pointed out that U n (z), etc., are one- valued functions of 
m — n + 1. The trouble is that the asymptotic expressions for U„(z) 
are many-valued functions of m even though U n (z) itself is not. 

Now that we have considered the saddle points t and U, we turn to 
a consideration of the paths of steepest descent in the /-plane which 
pass through them.* The path of steepest descent which passes through 
t , for example, is that branch of the curve 



Im[f(t)-f(to)] =0 



(10.8) 



for which k is the highest point (i.e., Re [f(t) - f(to)] ^ on it). The 



* Watson 22 has studied paths corresponding to Re(n) > when z is any com- 
plex number, and has given curves which are related to some of those shown in 
Section 11. 
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paths of steepest descent may be shown to have the following prop- 
erties : 

1. Let t = U + iU = r ex P (*0)- Then the paths of steepest descent 
either run out to t T = + °o with <,• — > Ira z or spiral in to t = as r = 
(constant) exp ( — m r 9/m,i). 

2. The steepest descent path through t may be computed by a graph- 
ical method based on* 



arg (dt) = arg t - arg (t - U) - arg (t - h). 



(10.9) 



,Vz2-2m 

ARG (z 2 -2ITl) 




Fig. 10.2 — Diiigriim showing the cuts in the complex m-plane, m — n + 1. 

If we draw the triangle to ti and bisect the interior angle at U by the 
line & fo then 



arg (dt) at to = angle tit o b . 



(10.10) 



If one goes clockwise in traveling from the side tot\ to tobo then arg c/£ is 
negative. Likewise, arg (dt) at h (on the path through h) is the angle 
between the side /i/ and the bisector /i&i of the interior angle at t\. 

3. When m has the critical value z 2 /2 the saddle points coincide: 
( = t y = z/2, and the paths of steepest descent start out from t = z/2 
in the three directions arg (t — z/2) = (arg z)/3 + 8 where 5 is 0, 2t/3, 
or — 27r/3. 

4. Some of the paths of steepest descent change their character as m 
goes from one region of the m-plane to another. This is illustrated in 
Section 11 for the case z = p exp (iV/4) where it is shown that the 

* A similar method was used in 1938 by A. Erdelyi in an unpublished study of 
the asymptotic behavior of confluent hypergeometric functions. 
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boundaries are given by 

Im\J{U) -Kk)\ =0, (10.11) 

or a similar equation involving another pair of saddle points, e.g., h and 
h exp (i2w). In this equation z is regarded as fixed and t , h are functions 
of m defined by (10.2). It should be noted that although (10.8) defines 
a path of steepest descent in the 2-plane, (10.11) defines curves (bound- 
aries of regions) in the m-plane. 

5. If m is such that the path of integration for a particular function, 
say U n (z), passes through both * and t h each one will contribute to 
the value of U n (z). Furthermore, if m is such that 

Re [/(*>) - f(k)] = 0, (10.12) 

to and h have the same height and the two contributions have a chance 
of cancelling each other and giving a value of zero for U n (z). Thus 
(10.12) or some similar equation defines the lines in the m-plane along 
which the zeros of U n (z), etc., (regarded as functions of m) are asympto- 
tically distributed. 

6. The lines in the m-plane defined by (10.11) and (10.12) may be 
obtained by substituting the values (10.2) for t and h in 

/(/o) - f(h) = k 2 - t{ - 2t n ti log (*,/&), (10.13) 

and setting the imaginary and real parts, respectively, to zero. How- 
ever, instead of dealing with m directly it is easier to use w = ii + iv 
defined by 

w = log (to/h) = log | to/h | + z'(arg * - arg h), (10.14) 

m = z 2 /(cosh w + i), (10.15) 

where (10.15) follows from (10.14) and (10.2). Then (10.13) becomes 
f(to) — f(fi) = m(sinh w — w) 

= 2 2 (sinh w - w) 0°- 16 ) 

cosh w + 1 

The inequalities (10.7) show that 

U ^ 0, | V | ^ 7T. 

7. For the special case z = p exp (iV/4), (10.16) gives 

(cosh u + cos v — v sin v) sinh u = (cosh u cos v + 1) w, (10.17) 
(cos v + cosh u + u sinh u) sin v = (cosh u cos v + 1) f> 
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respectively, for Im\f(t ) - fit,)} = and Re [/(&>) - /(*,)] = 0. These 
equations are plotted in Fig. 10.3. It will be noted that a curve is shown 
for v > t even though this puts w outside the allowed rectangle. This 
is done because one of the paths of integration, W, passes through both 
t and h exp ( — i2ir) when m is in a certain region, and the correspond- 
ing zeros of W n (z) lie on the curve defined by 

Re If (h) -f(hexp [-&*))] = 0. 
It may be shown that a curve corresponding to 
f(U) -fikexp |-t»r)) 

with — x < v < w may be obtained from the curve corresponding to 
f(to) — Kk) with ir < v < 3tt by simply subtracting 27r from v. This is 
done on Fig. 10.3. 




ARG m =-90° 



•0 1-90° 



BOUNDARIES 

*Im[f(t )-fCt,)] = o 

LINES OF ZEROS (EXCEPT CURVE (a)) 
'Re[f(t )-f(t,)] = o 



ARG m = 90° 



•O |90° 



/ Re[fao)-f(t,e-L2rr)] = o 



(a) 



ARG m = 270 



m — ■*•<> 270 



Fig. 10.3 — Boundaries of the regions shown in Fig. 11.2 and lines of zeros 
shown in Fig. 12.1 as they appear on the w = u + iv plane when z = ?' ,/2 p. 
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3. m = (0.54953 . . . )p ~ 0.55p 2 . Fig. 11.1 (b). This value of m 
marks the change in type of path. Since Im/(£i) = Im /(fo) is satisfied, 
the paths through /<, and h have the same equation [see (11.2)], and 
there is a chance for a situation like that at U in Fig. 11.1 (b) to occur. 
The high point of the path U is at fr, and it goes continually downhill 
on either side of h although its direction changes sharply by 90° at /<,. 
The point m = 0.55p 2 is just one point on the boundary between regions 
in the ra-plane corresponding to various types of paths. The boundary 
lines are obtained by solving condition (10.11) for m as outlined in 
Items G and 7 of Section 10. Mapping the boundary lines 

Im [/(*,) - f(hj\ = 

from the auxiliary to-plane (shown in Fig. 10.3) to the m-plane with the 
help of m = ip 2 /(cosh w + 1) gives the boundaries between the regions 
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I, II, III shown in Fig. 11.2. It may be verified that for large negative 
values of m,i the boundary lines in Fig. 11.2 are given approximately by 

m r = ± 2 3/ W | m | 1/2 . (11.4) 

There was a period, while these curves were being worked out, during 
which it appeared that the regions I, II, III told the entire story. How- 
ever, when small values of m were studied it was found that region I 
splits up into the two sub-regions, la and lb, such that the boundary 
between them is given by 

Im Wi) - Kke~ 2ri )\ = 0. (11.5) 



Vfi Wi. 




STEEPEST DESCENT /J. 

To^e- 1277 " 



Fig. 11.3 — Paths of steepest descent for | m \ = 0.005 p 2 , z = t'*' 2 p. Away from 
t = t l all paths look much like Fig. 11.3(a). 
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This is of the same form as (10.11). That h exp ( — 2id) is a saddle point 
follows from differentiation of the equation 

f(te~ 2ri ) = f(t) + 2wim r - 2irm,-. (11.6) 

Combining (11.5) and (11.6) shows that the boundary between la and 
lb is given by m r = 0. This is indicated on Fig. 11.2. 
We now examine the paths of steepest descent when m is small. Fig. 11.3 
(a) gives a large view of all the paths, irrespective of arg ?n, when m / p 
is small. 

4. m = 0.005p 2 . Fig. 11.5 shows the vicinity around t\. 

5. | m | = 0.005p 2 , arg m = ir/2 -0.05. Fig. 11.3 (b). 

6. | m | = O.OOSp 2 , arg m = t/2. Fig. 11.3 (b). After passing through 
t\ the path encircles the origin clockwise and runs down into the saddle 
point at t = t\ exp ( — 2iri). Since m» is positive, (11.6) shows that 
h exp ( — 2wi) is lower than fo. The path for arg m = r/2 — 0.05 sug- 
gests that from k exp ( — 2ti) the path runs out to °o exp (— iri) along 
the path of steepest descent which lies directly under (on the Riemann 
sheet for — St < arg t < — ir) the path which runs from h to 
t = co exp (i7r). It follows from (11.6) that, as t traces out a path of 
steepest descent through t\, t exp ( — 2-wi) traces out a path of steepest 
descent through t\ exp ( — 2-wi) directly under the path through t%. 

7. \m\ = 0.055p 2 , arg m = tt/2 + 0.05. Fig. 11.3 (b) shows that 
after passing through t\ the path of steepest descent spirals in to t = 0. 
According to (11.3), the spiral is given by 

r £H (constant) exp ( — m r 9/mi) (11-7) 

when r is small and 6 large. Two things are to be noted. First, the type 
of path is different from that for arg m = w/2 — 0.05. Hence arg m = 
tt/2 marks a change of type similar to that shown in Fig. 11.1 (b), except 
that here h exp ( — 2-wi) takes the place of t . Condition (11.5) takes the 
place of condition (10.11), and is satisfied by virtue of m r = when 
arg m = tt/2. 

The second thing to be noted is that up until now all of the paths of 
steepest descent have ended at ± °° and U, V, W could be deformed 
into them without difficulty. How can we deform U, for example, into 
a path of steepest descent when the path through ti spirals in to t = 0? 
The way to deal with this problem is shown in Fig. 11.4 where U is 
continuously deformed into two portions, one coinciding with the path 
through ti, as shown in Fig. 11.3 (b), and the other with the path of 
steepest descent through h exp (— 2ri). The second portion lies directly 
"underneath" the first portion. 
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In Fig. 11.4 the dashes mean, as before, that the path of steepest descent 
is on a sheet of the Riemann surface other than | arg 1 1 < t. The alter- 
nate dots and dashes are used to indicate that | arg t | > x and that in 
addition the path lies directly under the curve it parallels. Although in 
Fig. 11.4 the two kinds of dashed curves are joined at about arg t = 
— 3x — x/4, they actually should spiral in to t = before they connect. 

8. | m | = 0.005p 2 , arg m = x. Fig. 11.5. For arg m = x, (11.6) shows 
that U and h exp ( — 2xt) are of the same height. 

9. m = 0.005p 2 , arg m = 3tt/2. For x < arg m < Sir/2, k exp (-2xf) 
is higher than h and the paths spiral into t = counterclockwise. At 
arg m = 3x/2 the rate of spiralling is zero and we have the path shown 
in Fig. 11.5 (which is the path for arg m = x/2 rotated by 180 degrees). 
Here arg t\ = 5x/4. 

10. m = 0.005p 2 , arg m = —x/2. The paths for arg m equal to —x/2 
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Fig. 11.4 — Deformation of path of integration U into path of steepest descent 
through h when m = 0.005 p 2 exp (jV/2 + i0.05). 
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and 3ir/2 have the same shape and both are highest at the saddle point 
whose argument is — 3tt/4. In both cases the contributions are the same 
and hence the value of U„(z), for example, is the same for arg m = —ir/2 
as for Sir/2 (as it must be since our parabolic cylinder functions are one- 
valued functions of m). 

Before leaving the region around m = we point out that when 
| m/p | <JC 1 the path of steepest descent through / is almost inde- 
pendent of arg m. Also, the curves of steepest descent for 



i/r(w) - i / e l t 



dt 



(11.8) 



PATHS SUPERPOSED 
BUT ARGUMENTS OF t 
DIFFER BY 277 




'.STEEPEST 
> ASCENT 

H i 

Fig. 11.5 — Paths of steepest descent for \m\ = 0.005 p 2 , z = /" 2 p. These 
curves ure much the same as those in Fig. 11.3(b) except for the values of arg in. 
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behave in much the same way as those just described. The line m r = 
divides the w-plane into two regions corresponding to different types of 
paths, and the negative real axis is a line of zeros corresponding to 
Re Ufa) — f(ti exp (-2irt))] = where h — m is the saddle point. 

11. | m | = 0.55p 2 , arg m = w. Fig. 11.6 (a). This value of m marks a 
change in the type of path. 

12. | m | = p 2 , arg m = t. Fig. 11.6 (b). 

13. | m | = p\ arg m = tt/2, arg (ip 2 - 2m) = 3tt/2. Fig. 11.6 (c). 
The complication of the paths in Fig. 11.6 (c) is due to the superposi- 
tion of two boundaries in the ?n-plane. Im /(f ) = Im f(ti) accounts for 



Wfi 



U/* Vf.Wj. 










_i^— — 




a ^i^ 




m=/9 2 e i - 7r 


(b) ■' 







i 



ARG (Z 2 -2IT1) -3TT/2 





Fig. 11.6 — Paths of steepest descent for miscellaneous values of m with z = 
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the path running from k to k, and Im f(ti) = Im /[fo exp ( — 2ri)] for 
the one running from k to to exp ( — 2iri). The saddle points in order of 
their height are t , k, t exp (— 2*t), J being the highest. 

14. w = ip /2. Fig. 11.6 (d). Here k = k and the dashed lines go 
into the saddle point at k exp ( — 2iri). The paths of steepest descent 
change their directions upon passing through the saddle points. 

12. ASYMPTOTIC EXPRESSIONS FOR U„(z), V n (z), W„(z) 

The asymptotic expressions given here are for z = i p and z = f ' p 
[with i 1 ' 1 = exp (t7r/4)] when n is not too close to z 12. As mentioned 
earlier, there is a close relation between our results and those given by 
Schwid.* The main difference is that we regard n as variable and z 
as fixed while Schwid regards z as variable with n fixed. Another point 
of difference is that in place of the ra = n + 1 which appears in our 
expressions for t and k the quantity ft + M appears in Schwid's work. 
The quantity 2n + 1 appears to enter naturally when the asymptotic 
values are obtained from the differential equations. This is seen when 
the WKB method is applied to equation (9.7). 

By examining the paths of steepest descent shown in the figures of 
Section 11 we can determine the saddle points corresponding to U„(z), 
etc., (for z = i m p) for various values of ft. The contributions to the 
integral (10.1) from the saddle points / to h were discussed in Section 
10. The contribution from the saddle point k exp ( — 2iri) (which enters 
when z = t 1/2 p) is, from (11.6), exp (i2irm) times the contribution from 

k. 

Although we shall be concerned mainly with asymptotic expressions 
for the parabolic cylinder functions themselves, expressions for their 
derivatives may be readily obtained. Thus U' n ( z ) = dU„(z)/dz has the 
asymptotic expression 

U' n (z) ~ 2t [contribution of U to U„(z)] 

+ 2k [contribution of *i to U n (z)] (12.1) 

+ 2k [contribution of k exp ( — 2iri) to £/„(*)] 

and similar expressions hold for V' n (z), W' n (z). These follow when we 
note that differentiation of the integrals (9.1), which define the functions, 
introduces a factor 2t into the integrand. Of course, if the path of in- 
tegration does not pass throught a particular saddle point, its contri- 
bution to (12.1) is zero. Upon replacing * and k by their expressions 



Reference 20, page 478. 
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(10.2) and subtracting the corresponding expression for zU n (z) we 
obtain 

'UM = U' n (z) -z U n (z) 

~ {z - 2m) 1 ' 2 [(k contribution) - (t\ contribution) (12.2) 

— {t\e~ 2 * 1 contribution)] 

where 'U„{z) is the function defined by (4.19). The same is true for 
Tn(*) and 'W n (z). 

Consideration of the various paths of integration shown in Section 11 
leads to the results shown in Table 12.1. The leading terms of the 
asymptotic expansions are listed for the various regions of the ra-plane 

Table 12.1 — Leading Terms in the Asymptotic Expansions for 

U n (z), V„(Z), W„(Z) WHEN Z = i W p, p > 



Region in w-plane 
m = n + 1 


IM«'"V) 


F n (»"2p) 


Wnd^p) 


la 
II 
lb 
III 


A, 

Ai — A a 
(1 - i in )Ay 
(1 - ;<»)A, 


Ao 
Ao 
Ao 

A Q - Ai 


-Aq - Ai 
-Ax 
-Ao - Ai + i 4n Ai 
-A + z 4 »Ai 



shown in Fig. 11.2. If the next order terms are required, they may be 
obtained from (10.4) and (10.5). 
The notation used in Table 12.1 is as follows: 



•1/2 



z = i"~p, m = n + 1, i = exp (i-w/2), 
- 7T /2 < arg m ^ 3tt/2, - tt/4 < arg t Q ^ 3tt/4, 
- tt/2 < arg (ip 2 - 2m) ^ 3tt/2, - 3tt/4 < arg k ^ 5tt/4, 
to = [i m P + (ip - 2m) m }/2, h = [i m P - (ip - 2m) m }/2, 

A = [t ll2 (ip 2 - 2m)- l,i /2iT m } exp /ft), (12 . 3 ) 

A, = [h m {ip 2 - 2my ll4 /2T 112 ] expf(h), 



f(h) = zk + g ~~ m 



iog«, -|(i -log I -log I) + ,«■*, 

f(h) - rf, + f - m log (, = |(l - log I - log I) + i"W 
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Sometimes it is helpful to use 

/„ \-i/2 Y m ( -, i m 
(2tt) exp -(1 - log - 

fl/r(^^) = l/r(l 4- n/2) for 

—w/2 < arg m < tt/2, ( 12 - 4 ) 
r"T (L=^W = z-"- 1 r(-7i/2)/27r, 

tt/2 < arg m < St/2, 

where the last line is obtained by setting m exp ( — iri) for m in the 
second line. 

The asymptotic expansions for regions lb and III may be obtained 
from those for la and // by using equations (9.11) and (9.13). However, 
the work is more difficult than one might suspect at first glance. 

Incidentally, the leading terms in the asymptotic expansions (9.15) 
and (9.17), which hold when p — > *> and n remains fixed, may be ob- 
tained by considering the entries for la and lb in Table 12.1. 

It is sometimes convenient to use the limiting forms of the asymptotic 
expressions when | m | ^> p . In this case, for z = i p, 

2 



2'o -* z ± i(2mf 2 
2k->zT i(2m) m 



1 -2^I + 0(W " S ' !) - 



, - 2 -l+o (m -), 



(12.5) 



log h/t — > + Mr + iz(2m) 



-1/2 



2 + 



6m 



+ 0(m- 5/2 ), 



where the upper signs hold when — tt/2 < arg m < ir/2 and the lower 
ones when ir/2 < arg m < Sir/2. Substituting (12.5) in (12.3), neglect- 
ing the higher order terms, and setting 



r> r»-3/2 —1/2 

B = 2 7r exp 



[f(i-i*;)+^4 



a = exp [- p(2m/i) in ], 

«i = exp [p(2m//') 1/2 ] = V«o , 



(12.6) 



converts Table 12.1 into Table 12.2. 

In this table B, a Q , a h are defined by (12.6); m = n + 1; — 7r/2 
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Table 12.2 — Leading Terms in the Asymptotic Expansions for 

Un(z), V n (z), W n {z) WHEN Z = l' 2 p AND | 2n | » p 



Region in w»-plane 
(«« = n + 1) 


J/n(»" 2 p) 


V n (fii* P ) 


W,.(i V2 p) 


la 
II 

lb 

III 


i n Bu\ 

(1 - i*»)i-»Ba 

(1 - {^{-"Bao 


—i- n Ba 

— i~"Bao 

i n B ai 

B(i n ai — i~ n ao) 


B(i~"oc — i n ai) 

— i"Bai 
B(i 3n ao — i~"ao 
— i n ai) 
B(i 3n a — i"oti) 



< arg m ^ 3ir/2; and in regions la and lb arg m is approximately 
— 7r/2 and Sir/2, respectively. Gamma functions may be introduced 
into the expression for B with the help of (12.4). It may be verified that 
the functions do not change, except for negligible terms, in crossing 
over the boundary from la to 76 (a and ai are interchanged and B is 
changed by the factor exp (— m-wi)). 

Since the zeros of our functions, regarded as functions of n, occur 
(asymptotically) when the contributions from two saddle points cancel 
each other, we may look at Table 12.1 and pick out regions which may 
possibly contain zeros. Thus, A may equal Ai along the line | A \ = 
| ili |, i-e. very nearly Re[f(t ) — f(ti)] = 0, in the w-plane. These lines 
were discussed in Item 7 of Section 10 and are plotted on the auxiliary w- 
plane in Fig. 10.3 When plotted on the m-plane the lines appear as 
shown in Fig. 12.1 The condition Re[f(io) — f(h exp ( — 2iri))] = gives 
the line I Aq 1 & I t 4n Ai I for some of the zeros of W„(i l2 p). 



y* (z) 




ZEROS OF U n (Z) 
I = l/J 2 /2 



Un (Z) 



|W n (z) 



.-^-~ W n (z) 

Fig. 12.1 — When U„(z), V„(z), and W n (z) are regarded as functions of n their 
zeros lie on the lines indicated when z = i ll2 p. The three branches coming out from 
in = ip 2 /2 are lines along which | Ao | = I A\ \ and the branch for W„(z) coming 
down from m = is a line along which I Ao j ~ 
in Table 12.1. 



i in A\ | where Ao and Ai appear 
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The location of the zeros far out on the lines of Fig. 12.1 may be ob- 
tained by writing the appropriate expressions of Table 12.2 as propor- 
tional to B times a cosine or sine. Examination of the trigonometrical 
terms shows that 

U n (i ll2 p) has zeros at n « 2/c + 1+ i'%k m /ir, 

V n (i m P ) has zeros at n « -2k + i 3/2 4pfc 1/2 /ir, (12.7) 

W n (i ll2 P ) has zeros at n « -2k + r 1/2 4p/c 1/2 /T, 

where k is a large positive integer. Of course, U„(z) also is zero when n 
is a negative integer. 

So far we have been dealing with z = % l2 p. Now we consider the case 

—1/2 
Z = I p. 

Asymptotic expressions which hold when z = i~ m p may be obtained 
from Table 12.1 by using the relations (9.11) between functions of z 
and of its complex conjugate z*. Thus, for example, V a +ib(i~ m p) is equal 
to the complex conjugate of W a - ib (i m p). These relations, and relations 
such as 

[t for z = i m p, n = a - ib]* = to for z = r m p, n = a + ib 

[f(to) for z = i m P , n = a - ib}* = j{h) for z = T 1/2 p, n = a 4- ib 

(12.8) 

have been used in constructing Tables 12.3 and 12.4 from Tables 12.1 
and 12.2 The interchange of V n (z) and W n {z) should be noted. The 

TTl- PLANE 

m = n+i 




W n (z) 



m' 



m = -L/> 2 /2 

.ZEROS OF Un (Z) 



n' 



Fig. 12.2 — Regions in the complex m-plane corresponding to different asymp- 
totic expressions when z = i~ U2 p. The lines on which the zeros of the various 
functions lie are shown by the dashed lines. The corresponding information for 
2 = i ll2 p is shown on Figs. 11.2 and 12.1. 
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Table 12.3 — Leading Terms in the Asymptotic Expansions for 

U n {z), V n (z), W n (z) WHEN Z = r ll2 p, p > 



Region in »t-plane 
m = n + 1 


U n (i-llt P ) 


K„(*-i'*p) 


W n (i- U, p) 


la' 
II' 
lb' 
III' 


t A[ 

Ai Aa 

(1 - i-*»)Al 
(1 - i- in )Al 


An Ai 
-Ai 
-Ai - (1 - i-*)Ai 
-A» + i-*»A[ 


A'o 

A' 

, A 'o , 
Ao A\ 



regions in the m = n -\- 1 plane corresponding to the different asymptotic 
expressions are shown in Fig. 12.2. The boundaries are simply those of 
Fig. 11.2 reflected in the real ra-axis. The lines of zeros are also shown in 
Fig. 12.2, and are reflections of those of Fig. 12.1 except for the inter- 
change of V n (z) and W n (z). 

Table 12.3 may also be constructed by returning to the paths of in- 
tegration shown in Section 11. It may be shown that corresponding to 
every path of steepest descent for z = i 1 2 p, n = ri\ there is another 
path, obtained from the first by reflection in the real f-axis, which gives 
the path of steepest descent for z = i~ 2 p, n = n*. 

The notation used in Table 12.3 is as follows: 



z = i 



-1/2 



'p, m = n + 1, i — exp (iV/2), 
-3tt/2 g arg ra < tt/2, -3tt/4 ^ arg k < ir/4, 
-3tt/2 ^ arg (-ip 2 - 2m) < tt/2, -5tt/4 ^ arg k < 3tt/4, 
k = \T m P + {-ip 2 - 2m) 1/2 J/2, k = [r U2 P - (-ip 2 - 2m) ll2 ]/2, 



Ao = [k ll \-ip - 2m)- 1,i /(-2iw 1 ' 2 )] exp /ft), 
A[ = [ti'\-ip* - 2m)- ll4 /2T ll2 ]expf(t 1 ), 



(12.9) 



m 



m 



m 



to 



f(to) =zt + '-^-m log h = S 1 - ^g £ - log f + C** , 



.-1/2. 



m 



to 



f(ti) = zh + - - m log k = g-l 1 - log - - log i- ) 4- -T^pti . 
Sometimes it is helpful to use 



k 



-—1/2. 



(2tt)" 1/2 exp 



;"(i-. 0.=)] 



Vr(^±i 



= l/r(l + n/2) for -tt/2 < arg m < tt/2, (12.10) 



tT f 1 m j / 2t = i n+1 r(-n/2)/27r, -37r/2 < arg m < - tt/2. 
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Table 12.4 — Leading Terms in the Asymptotic Expansions for 

Un(z), V n (z), W n {z) WHEN Z = C l '*p AND | 2» | » p 2 



Region in m-plane 
m = n + 1 


t/n(«-"V) 


p»(t-" J p) 


IP„(*-i'»p) 


la' 
II' 

lb' 

III' 


i-»B'a[ 
B'(i n ao + i~"ot'i) 
(1 - i-*")i n B'a' 

(1 - i-*")i"B'ao 


B'(i n ao — i~ n a[) 

-t-»B'«! 

B'(i~ 3n a' — l B ato 
B'(i- 3n o:o - *-»a0 


— {"B'a'o 
-i"B'a' 
i~ n B'ai 

B'(i~ n ct[ — i n a' a ) 



The notation used in Table 12.4 is as follows: 

m = n + 1, i = exp (ttr/2), -37r/2 ^ arg m < tt/2, 

B ' = 2-V-exp^(l-.o^)-i-], 

«o = exp[-p(2im) 1/2 ], (12.11) 

a[ = exp [p(2iwi) 1/2 ] = l/«o, 

B' may be expressed in terms of gamma functions with the help of (12.10.). 
Approximate expressions for the zeros are given by the complex 
conjugates of (12.7). For example, if k be a large positive integer such 
that 2k » p 2 , the zeros of W n (f 1,2 p) are at n = n(k) where ■T n a[ = exp 
(iirk) and 

n(k) fist -2k + r 3/2 4pfc 1/2 A - 4ip 2 /7r 2 . (12.12) 

Here the approximation has been carried out one step further than in 
(12.7) We also have for the quantities in (7.11) 

[dW n (r ll2 p)/dn) n=n{k) « (-) k+1 B'i[* - p(i/k) m ), (12.13) 

[* 2 'T n (r 1/2 p)/sin irn) n=n(k) « (-) k+1 2B'i. 

13. ASYMPTOTIC EXPRESSIONS FOR U„(z), ETC., WHEN W IS NEAR Z 2 /2 

The asymptotic expressions given in Section 12 fail when n is near 
2 2 /2. Expressions for the parabolic cylinder functions which hold for 
this region have been given by Schwid.* More recent studies of this sort, 
based on differential equations, have been made by T. M. Cherry 4 and 
F. Tricomi. 25 Their results suggest the possibility that our expressions 



* Reference 20, page 478. 

24 Uniform Asymptotic Expansions, J. Lond. Math. Soc, 24, pp. 121-130, 1949. 
Uniform Asymptotic Formulae for Functions with Transition Points, Am. Math. 
Soc. Trans.,' 68, pp. 224-257, 1950. 

16 Equazioni Differenziali, Einaudi, Torino, pp. 301-308, 1948. 
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for the electromagnetic field which contain Airy integrals may be re- 
placed by more accurate, but also more complicated, expressions. In 
dealing with our functions we shall work with the integrals and our 
procedure is somewhat similar to that used by Rydbeck.* First however, 
we point out that when we write (as suggested by the work of Cherry 
and Tricomi) 

y = e~' 2 ' 2 T n (z), ax = z-(2n+ 1) 1/2 , 2(2n + l) 1/2 a 3 = 1, (13.1) 

the differential equation (9.7) for the parabolic cylinder functions goes 
into 

Q-zy- 2- i,3 (2n + D^W (13.2) 



dx 

The Airy integrals Ai(x) and Bi(x) (and also Ai[x exp (± i2ir/2>)]) 
discussed later in this section are solutions of 

g - „ - 0, (13.3) 

and therefore we expect that approximate solutions of (13.2) are given 
by, for example, 

y = aA.i(xT)[l + 0(rT 2/3 )] (13.4) 

where the 0(n -2/3 ) term corresponds to the particular integral of (13.2) 
when the y on the right hand side is replaced by its approximate value 
Ai(xi"). Here Ci is independent of x (or z) but may depend on n, and v 
may be or ±4/3. 

Since the labor of computing C\ is considerable, we shall work out the 
approximations directly from the integrals. 

We shall consider the case z = i ll2 p, p > 0, first. When n + 1 = to 
= mo = ip 2 /2 the saddle points to and h coincide at t 2 = i 1,2 p/2. Con- 
sequently only those portions of the paths of steepest descent which 
lie near U are of importance. This is true even if m is not exactly equal 
to mo. We therefore regard 

f(t) = -t 2 + 2zt - m log t (13.5) 

in (10.1) as a function of the two variables t and m (linear in m) with z 
fixed at i ll2 p. Expanding (13.5) about t = ti , m = mo gives 

r m Wo , Wo (w — Wo) , TO 
T - T log T -^^log T _ (i36) 

- 4:(t - hf/Zz - 2(w - m )(t - k)/z + • • • 



m = | + 



* Page 87 of Reference 21 cited on page 478. 
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where we have used 

t 2 = 2 /2 = i m P /2 = (m /2) m (13.7) 

and have arranged the terms within the brackets so that they represent 
the first two terms in the expansion of 



S-2l<«?~l«*[<3*)"/- /m + 1 



(13.8) 



2 2 °2 

about wio . 

The paths of steepest descent in the i-plane when m = m are shown 
in Fig. 11.6(d). The three branches start out from t = t 2 in the directions 
arg (t — t 2 ) = 15°, 135°, and — 105°. In this section we take the paths of 
integration to be those of Fig. 11.6(d) even when m is not exactly 
equal to m . Since we are dealing with asymptotic expressions we may 
confine our attention to the region around t = t 2 where the paths of 
integration are essentially straight lines [the contributions from t 2 exp 
(— 2iri) are negligible]. 

When (13.6) is set in the integral 

(1/2W) / exp [/(/)] dt (13.9) 

we see that the initial directions of the branches are such as to make 
(/ — k) 3 /z positive (arg z = 45°). Some study of (13.6) and of the Airy 
integrals we intend to use suggests that we change the variable of inte- 
gration from i to s and introduce the parameter b where 

l -k = s(z/4) I/3 , b = (m - m„)(2/2 2 ) 1/3 = (m - m )/m m . (13.10) 

This and (13.6) converts the integral (9.1) for V n (i m p) into 

z(27r) 1 ' 2 r^— — J ( 13-11) 

exp [ — bs — s/3 + • • •] ds. 
When we use the Airy integral defined by 



Ufa) = T ~ l [ cos fat + t 3 /S) dt 

Jo 

= (i m /2w) [ exp [—T m xs - s 3 /3] ds, 

Joo exp (i2x/3) 



(13.12) 



cp («t/3) 

we obtain 

vl/3 «Vl 






(13.13) 
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In order to obtain expressions corresponding to (13.13) for c/„(z), W„(z) 
we examine Fig. 11.6(d). We have already seen that the limits of in- 
tegration for s, in the integral (13.11) for V„(i m p), are 
[oo exp (i2ir/3), oo]. In the same way it follows that the limits for 
U n (i in P ) and W n (i m p) are [« exp (- z'2tt/3), «, exp (*2ir/3)] and 
[oo , oo exp (— i2r/3)] , respectively. When we take s' = s exp (+ i2t/3) 
as new variables of integration (with the upper sign for C/„(z) and the 
lower one for W„ (z)), the integrals corresponding to (13.11) go into 
Airy integrals. 

We can write our results for z = i p, when n is close to ip /2, as 
follows: 

U n (i}' 2 P ) ~ Ci^AiQn" 9 ), 

V n (? ,2 p) ~ Ci- 7l6 Ai(bi 213 ), (13.14) 

W n {i l2 p) ~ Ci m AiQn- 213 ), 



where 



\l/3/ \'/2„«P 2 /2/ T ,/ m + 1 



C = (p/4) lM (2r) ,/ V' , ' /1 /r i ., 

& = (2/p 2 ) 1/3 (m - ^p 2 /2)^- 1/ ^ 
i = exp (it/2), m = n + 1. 



(13.15) 



The asymptotic expansions whose leading terms are given by (13.14) 
may be obtained by the method used by F. W. J. Olver 26 to study 
Bessel functions. 

Ai{x) and its derivative have been tabulated for positive and negative 
values of x* Here we shall use the definitions and results as set forth in 
Reference 11. These tables and (13.14) enable us to obtain values of 
Un{i m p) along the rays in the ?n-plane defined by arg (m — ip /2) = 7r/6 
and — 5ir/Q. Along the x/6 ray bi 6 ' 3 is negative. Since the tables show 
that the zeros of Ai(x) occur when x is negative, it follows that the 
zeros of C/„(i 1/2 p) occur on the x/6 ray. In the same way it is seen that the 
zeros of V n {i w p) and W n {i m p) occur on the 5ir/6 and the — ir/2 rays, 
respectively. This agrees with Fig. 12.1. 

The Airy integral defined by 



* Reference 11, page 424. 

26 Some New Asymptotic Expansions for Bessel Functions of Large Orders, 
Proc. Cambridge Phil. Soc, 48, pp. 414-427, 1952. 
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Bi(x) = 



—2/3 p r »exp (»2»/3) 

V 



2*7T |_Joo exp (— *2jt/3) 



£ 



+ 

.«c P (-»2t/3)J 



.—2/3 3 



exp \-r l "x8 - a*/Z] ds (13.16) 



7T Jo 



t* A , . /t 3 



exp ( -- + xtj 4- sin I - + art 



d/ 



is also tabulated in Reference 11 where it is shown that 

Ai(xi 41 *) = i w [Ai(x) - iBi(x)]/2, 

(13 17) 
Ai{x€ m ) = C*"[Ai(x) + iBi(x)]/2. " ' 

With the help of these relations we may evaluate the expressions (13.14) 
for U n (i 2 p), etc., on any one of the six rays 

arg (m - ip 2 /2) = ±5tt/6, ±/tt/2, ±tt/6. 

When b is a general complex number the expressions (13.14) may be 
evaluated with the help of the modified Hankel functions hi(a), h 2 (a) 
tabulated in Reference 27 for complex values of a. The relation needed 

is 

Ai(a) = -L l h ( -a) + -L M -«), 

2A; 2/c* ( 13#18 ) 

k = (12) 1/6 2- 1/3 . 

When | arg a \ < ir we have the asymptotic expansion 

Ai(a) ~ 2~V _1/ V 1/4 (exp [- (2/3)a 3/2 ])(l - 5/48a 3/2 + • • •) (13.19) 

and when | arg (— a) | < 27r/3 we have 

Ai(a) ~ 7T~ 1/2 (- «) _1/4 sin [(2/3)(- a) 3/2 + tt/4] . (13.20) 

Both of these expansions follow from the discussion of the asymptotic 
behavior of hi(a) and /12(a) given by W. H. Furry and H. A. Arnold. 27 

Asymptotic expressions for U n (i~ l 2 p), • • • valid when n is near —ip/2 
may be obtained by applying the relations U„(z*) = [U n *(z)]* • ■ • given 
by (9.11) to the expressions (13.14) for U n (i 2 p), • • • : 

£/ n (r 1/2 p) ~ C'i~ m Ai(b'i~ m ), 

F.(f l/i P ) ~ C'r m At(b'r' 3 ), (13.21) 

TF n (r 1/2 p) ~ C'i V6 Ai(b'i~ 213 ), 

27 Tables of the Modified Hankel Functions of Order One-Third and of Their 
Derivatives, Harvard Univ. Press, 1945. 
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where 

c = (p/^'WV^Vr (*±!) , 

V = (2/p 2 ) 1/3 (m + ip 2 /2)i 1/3 , (13 ' 22) 

i = exp (iir/2), m — n + 1. 

In (13.14) 6i 6/3 = -6 and in (13.21) &V 8 ' 3 = -6' since 4*'(a) is a 
single- valued function of a. It is interesting to note that the factor 
i 116 in the expression for U n (i m p) gives the direction of that one of the 
three paths of steepest descent (in the i-plane) which is not traversed 
in getting U n (i ll2 p). The same sort of thing is true for the remaining 
expressions in (13.14) and (13.21). 

The functions 

'U n (z) = exp (z 2 /2)d[U n (z) exp (- z 2 /2)]/dz, 

denned by (4.19), may be computed from (13.21) when z = f~ m p. We 
need the relations d/dz = i 2 d/dp and 

m = -ip 2 /2 + b'(p 2 /2i) m , 

(13.23) 
db'/dp = (2/3)(2* P ) 1/3 (; - m/p 2 ) = i(2i P ) 113 - 2&'/3p, 

which follow from the definition of 6'. When the differentiations are 
carried out we obtain 

'U n {r m p) ~ (2p) U3 C'i- ll3 Ai , (b'i f - 6,z ), 

'V n (i- ll2 p) ~ {2p) ll3 C'i w Ai'{Vi 2lz ), (13.24) 

'W n {r m P ) ~ (2p) 1,3 C'i 7,3 Ai'(Vi- 2 ' 3 ), 

In these expressions the prime on the Airy integral denotes its deriva- 
tive: 

Ai'(a) = dAi(a)/doc. (13.25) 



